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33 r Ccyih 

Enter the value of the mass (kg) : 1 
m = 

1 

Enter the value of the darrping coefficient (N.s/m) : 4 
c = 

4 

Enter the value of the spring constant (N/m) : 1 
k = 

1 

Enter the value of the initial position (m) : 0 
xl = 

0 

Enter the value of the initial velocity (m/s) : 1 
yl = 

1 


ans = 

Natural period 

T = 


6.2832 
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2~3 f Oc.-if 


Enter the value of time incrsnent in seconds ( < T as above ) : .5 
dt = 


0.5000 


Enter the value of the initial time in seconds (0) : 0 
tl = 


0 


Enter the value of the final time in seconds : 10 
tf = 


10 


ans = 

Time Displ. Vel. 


ans = 


0 

0 

1.0000 

0.5000 

0.1667 

0.2526 

1.0000 

0.1953 

0.0360 

1.5000 

0.1856 

-0.0235 

2.0000 

0.1667 

-0.0368 

2.5000 

0.1471 

-0.0371 

3.0000 

0.1290 

-0.0339 

3.5000 

0.1130 

-0.0301 

4.0000 

0.0988 

-0.0264 

4.5000 

0.0864 

-0.0231 

5.0000 

0.0756 

-0.0203 

5.5000 

0.0661 

-0.0177 

6.0000 

0.0578 

-0.0155 

6.5000 

0.0506 

-0.0136 

7.0000 

0.0442 

-0.0119 
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Displacement 


A - 3 f 



7.5000 

0.0387 

-0.0104 

8.0000 

0.0338 

-0.0091 

8.5000 

0.0296, 

-0.0079 

9.0000 

0.0259 

-0.0069 

9.5000 

0.0226 

-0.0061 

10.0000 

0.0198 

-0.0053 


ans = 


maximum arrplitude 


A = 


0.1953 


0.2 

0.18 

0.16 

0.14 

0.12 

0.1 


Displacement Vs. Time 


o 08 





» ^ .? > Ccn ~f~ 

Enter the value of the ness (kg) : 1 

m = 

1 

Enter the value of the danping coefficient (N.s/m) : 
c = 

1 

Enter the value of the spring constant (N/m) : 1 
k = 

1 

Enter the value of the initial position (m) : 0 
xl = 

0 

Enter the value of the initial velocity (m/s) : 1 

yl = 

1 


ans = 

Natural period 

T = 

6.2832 
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oL-Jt Cont 

Enter the value of time incranent in seco nds ( < T as above ) 
t = 

0.5000 

Enter the value of the initial time in seconds (0) : 0 
1 = 

0 

Enter the value of the final time in seconds - 20 
;f = 

20 


Displacemen Vs. Time 




^ -.5 f- Co a f 
Time Displ. 


0 

0.5000 

1.0000 

1.5000 
2.0000 

2.5000 
3.0000 

3.5000 
4.0000 

4.5000 
5.0000 

5.5000 
6.0000 

6.5000 
7.0000 

7.5000 
8.0000 

8.5000 
9.0000 

9.5000 
10.0000 

10.5000 
11.0000 

11.5000 
12.0000 

12.5000 
13.0000 

13.5000 
14.0000 

14.5000 
15.0000 

15.5000 
16.0000 

16.5000 
17.0000 

17.5000 
18.0000 

18.5000 
19.0000 

19.5000 
20.0000 


0 

0.3776 

0.5340 

0.5259 

0.4196 

0.2743 

0.1332 

0.0219 

-0.0499 

-0.0838 

-0.0882 

-0.0739 

-0.0510 

-0.0272 

-0.0076 

0.0058 

0.0128 

0.0146 

0.0129 

0.0093 

0.0054 

0.0020 

-0.0005 

-0.0019 

-0.0024 

- 0.0022 

-0.0017 

- 0.0010 

-0.0004 

-0.0000 

0.0003 

0.0004 

0.0004 

0.0003 

0.0002 

0.0001 

0.0000 

- 0.0000 

- 0.0001 

- 0.0001 

- 0.0001 


1.0000 

0.5182 

0.1260 

-0.1363 

-0.2692 

-0.2980 

-0.2580 

-0.1840 

-0.1036 

-0.0349 

0.0136 

0.0404 

0.0488 

0.0446 

0.0334 

0.0202 

0.0082 

-0.0006 

-0.0058 

-0.0079 

-0.0076 

-0.0060 

-0.0038 

-0.0018 

- 0.0002 

0.0008 

0.0012 

0.0013 

0.0011 

0.0007 

0.0004 

0.0001 

- 0.0001 

- 0.0002 

- 0.0002 

- 0.0002 

- 0.0001 

- 0.0001 

-0.0000 

0.0000 

0.0000 


Vel. 
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t ' J f' Oosi f~ 


maximum anplitude 

A = 


0.5340 

» (~c> b (Q ) 

Enter the value of the mass (kg) : 1 

m = 

1 


Enter the value of the danping coefficient (N.s/m) : 2 
c = 


2 


Enter the value of the spring constant (N/m) : 1 
k = 


1 


Enter the value of the initial position (m) : 0 
xl = 


0 


Enter the value of the initial velocity (m/s) : 1 
yl = 


1 


ans = 


26 



x. ~ 3 t- 0 c '^-r- 


Natural period 

T = 


6.2832 


Enter the value of time increment in seconds ( < T as above ) : .5 
dt = 


0.5000 


Enter the value of the initial time in seconds (0) : 0 
tl = 


0 


Enter the value of the final time in seco nds : 20 
tf = 


20 

ans = 

Time Displ. Vel. 

ans = 


0 

0 

1.0000 

0.5000 

0.3021 

0.3047 

1.0000 

0.3666 

0.0016 

1.5000 

0.3337 

-0.1103 

2.0000 

0.2699 

-0.1344 

2.5000 

0.2047 

-0.1225 

3.0000 

0.1491 

-0.0992 

3.5000 

0.1055 

-0.0752 
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cS - J 7- oca f- 


4.0000 

0.0732 

-0.0548 

4.5000 

0.0500 

-0.0388 

5.0000 

0.0337 . 

-0.0269 

5.5000 

0.0225 

-0.0184 

6.0000 

0.0149 

-0.0124 

6.5000 

0.0098 

-0.0083 

7.0000 

0.0064 

-0.0055 

7.5000 

0.0042 

-0.0036 

8.0000 

0.0027 

-0.0024 

8.5000 

0.0017 

-0.0015 

9.0000 

0.0011 

-0.0010 

9.5000 

0.0007 

-0.0006 

10.0000 

0.0005 

-0.0004 

10.5000 

0.0003 

-0.0003 

11.0000 

0.0002 

-0.0002 

11.5000 

0.0001 

-0.0001 

12.0000 

0.0001 

-0.0001 

12.5000 

0.0000 

- 0.0000 

13.0000 

0.0000 

- 0.0000 

13.5000 

0.0000 

- 0.0000 

14.0000 

0.0000 

- 0.0000 

14.5000 

0.0000 

- 0.0000 

15.0000 

0.0000 

- 0.0000 

15.5000 

0.0000 

- 0.0000 

16.0000 

0.0000 

- 0.0000 

16.5000 

0.0000 

- 0.0000 

17.0000 

0.0000 

- 0.0000 

17.5000 

0.0000 

- 0.0000 

18.0000 

0.0000 

- 0.0000 

18.5000 

0.0000 

- 0.0000 

19.0000 

0.0000 

- 0.0000 

19.5000 

0.0000 

- 0.0000 

20.0000 

0.0000 

- 0.0000 


ans = 

maximum arrplitude 


0.3666 
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Displacement Vs. Time 
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ln(X) 
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+ o +*'('+Jtf]* * — c *'* 

(*•,++ ZA*)* +[4, + 4fi+af J x + za 

5 - z f\^. ce - ^ ya^lf ] 7*^ >■ 
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- O 
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“ ^Ae 'mt/a./ Condi i/orts %/°) - a.noL 

£(o) — O j so/ut/on /»us/ b co/ts/dered f°r <s<a - c ^ 

^ l^uri/tq //e fc^/ Aa./f cyc/e^ u///A fn. 

Ah at 


are.. 


half Cyde. ~Dur,/,j // € fr?/ Aa./f cy c/e , -** 

moving from r><j U /o left, x(/) cf ru?y <zt/e So 
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%/t) - A sfn* m t + 3 Cosoj m t •+ 
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%Ct) = CSA^cJjr -tZ> Cos coj: — 

~ — 4J 2>SEh. 

If/e ayain meerure time t = 0 d//be bey nnlry of/be 2 half 
Cyc/e, w/A neuJ initial Condi//ortS found from Aeywed/ous 
equations <iS %(°) - (- * D +- Z/ J^) ancL ^ (e> / ~ ° 

7) - -% -t 3 mF and C - O 

^ R 

'.There* acfua/iorxfor /he second Aa/f Cye/e 

*(t) = (~/o + 3 /£) Ccns (*„i - f*£ 

%(t) - ~ (-Ze T 3 H* t 
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1. {l 



EfUah'o* of pfo fr /or] 

MX- -kepX - CX s=.> mx + CX t /c e jX=0 

^ ri j -^=- -!_ Cset proUtiH 2.SO) fs tit 

f A, t l/h z 

ttttcbut. Y r/ ^/ £o/;sfa/r(r for the. sysle/v. 
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-3-4 


Square Ey, (3.I-7) 




(JL t 

/ 

\x 0 J 

(/-/L 

(JL) X s 

9.U-A,) - 45 ** 

\ xj 

(denmmj x 


y. = >y t 

A =(“/».)* 
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J P 


ftU* 
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3^.J./-/7 /V^j = 
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■2.<*/fe 


[,-A* J* + [i ? A]‘ 

R + t I 


■^Ttiis is the function file for the imaginary port of prob. 3-6 

function [iiragp]=iirp(r) 

zeta=0.01; 

(2*zeta. *r). / ((1-r. *r).~2+ (2*zeta. *r) .~2) ; 
and 


%This is the function file for the real part of prob. 3-6 

function [realp] =rp(r) 

zeta=.01; 

realp=(l-r.*r). / ((l-r.*r) . /x 2+(2*zeta.*r) ,~2) ; 
ord 
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To boJancf. y rotate tAe sec.tor- * - d.OSf , /2J,C c/*e^«w« 

And mu/t/pty bp * Tie poft/o* of fie la/anC"?J ,s 

found Ay roU + 'mj TU Mo./ ve,jM /ZJ.f e.W. Mom 
pe$'t' 0 » (of ‘fl.S'-Mont r^f.) and redue/nj >/s ^4<f /5«« 
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3-IQ 



OrijinaJ unbcdance cw.f’rom Ref'. 

WlM tr/'oL> k/e/fJJ W t ~2 @ J35* cw " » *#*> reo.d,n^ 4.50 <g /#"*,. 

Resultant =. -Atfi" <£> 7(,S° fn** hereon.to./ 

To b*.Jance. rcUte Wj- 76 .£* And reduce to * - A3f oz. 

** * ^* 5 " 

Student? sbouJd c.Aeck wWA pnob. S~~^ 
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3-/2 The method used to balance the single disk, can 

be modified to balance the long rotor ahown In Pig P3-12 (a). 


Pig P3-12 (a) 

ntartfUM fttrtndF 



(1) First measure the original unbalance vectors w.r.t. a fixed 
reference on the rotor at the near and far ends as N and P 
respectively. 

(2) Next, a trial weight is added to the near end and the 
new vectors N 1 and P 1 , shown in Pig P3-12 (b) are recorded 

with W t1 . The dotted vectors A and o(A are due to W tl alone. 



Pig P3-12 (b) 


and since they are proportional to W tl , they can be related to 
e&tA other by a complex influence coefficient, or operator o< . 

(3) Next, the trial weight W t1 is removed and a new trial 
weight W ^2 is added to the far end. Measurements for Ng 
and Pg are then recorded, as shown in Pig P3-12 (c). 




Pig P 3-12 (c) 
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and 


Izll CjszL 

Again, the dotted vectors are proportional to W t2 • 
they are related to each other by another complex operator ^ 
(4) As in the single dish problem, the two dotted vectors for 
the near end must be rotated and combined to cancel the 
original unbalance N. Also the final balance weights 
and W b2 must be related to the trial weights W t1 and W t2 
by two additional complex operators Q and (fi 

*bi = e » t1 

w b2 = f W t2 

Thus to balance the long rotor, the following two equations 
must be satisfied 

9 A + <$>(?&) = - V 
0B ■+ 9 i^A) ~ "F 


Solving these for 9 and 


* 


e ® 



./I F =JL- 
C /- A 

otN — F 




we obtain 


In these equations, N and F are the original unbalance 

vectors, and and are obtained by the division 

oc and A = && . With 6 and 0 known, the 

A ' B r 

balancing weights W fe1 and W b2 are found. 
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Eq.3 .i-& 


'Xtt) = 


4J — T°it „ , _ ^eui~ 

x,e s^Ozrvf-ty+Tfgr.Z3r7^r ' 

•tr&nf/crit Sfe.«.*<.v 
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22-141 50 SHEETS 

22-142 100 SHEETS 

22-144 200 SHEETS 
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3 -S<e 


clear 

global m k a F w 
'Problem 3-5 (f. Solution' 
x0(1 )=0;x0(2)=0;t0=0; 

m=input('Enter the value of the mass (kg) : ') 

k=input('Enter the value of the spring stiffness : ') 

a=input('Enter the value of the damping coefficient a : ') 

F=input('Enter the amplitude of the forcing function : ') 

nc=input('Enter the desired number of cycles: ') 

w1=input('Enter the lower bound of the forcing frequency w : ') 

w2=input('Enter the upper bound of the forcing frequency w : ') 

iw=input('Enter the desired increment of the forcing frequency w : ') 

'Wait .... this might take some time' 

t0=0; 1=0; 

for om=w1:iw:w2 

tf=nc*2*pi/om; 

1=1+1 ;wf(l)=om;w=om; 
clear t x ampl 

[t,x]=ode45('p350d',t0,tf,x0); 

mx=size(x);n=mx(1,1); 

sum=0;j=0; 

for i=1:n-1 

if (x(i+1,2)<0 I x(i+1,2)==0) & x(i,2)>0 

j=j+i; 

ampi(j)=x(i,1 )+(x(i+1,1)-x(i,1))/(x(i+1,2)-x(i,2))*(-x(i,2)); 
sum=sum+ampl(j); 
end 

end 

avamp(l)=sum/j; 

end 

plot(wf,avamp);xlabel('Frequency w');ylabel('Amplitude average') 
title('Problem 3-5£- solution') 


function xdot=p35£d(t,x) 
global m k a F w 
if ,x(2)>0 I x(2)==0 
xdot=[x(2); 

(1/m)*(-a*x(2) A 2-k*x(1)+F*sin(w*t))]; 

else 

xdot=[x(2); 

(1/m)*(a*x(2) A 2-k*x(1)+F*sin(w*t))]; 

end 
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'Problem 3-57 Solution' 
clear 

xO=input('Enter the value of the initial position of the ball >0 : ') 
uO=input('Enter the value of the initial velocity of the ball : ') 
N=input('Enter the desired number of hits : ') 
w=input('Enter the value of the frequency of the table : ') 

'Wait...' 

inc=(2*pi)/(300*w); 

g=9.81; 

t0=0;t(1 )=0;x(1 )=x0;u(1 )=u0;y(1 )=0; 

j=i; 

for i=1:N 

while x(j)>y(j) I x(j)==y(j) 

j=j+1; 

t(j)=t(j-1)+inc; 

x(j)=x0+u0*(t(j)-t0)-0.5*g*(t(j)-t0) A 2; 

y(j)=sin(w*t(j)); 

u (j)=u0-g*(t(j)-t0); 

end 

t(j)=t(j-1 ) + (y(j-1 )-x(j-1 ))*(t(j)-t(j-1 ))/((x(j)-y(j))-(x(j-1 )-y(j-1))) 
thit(i)=t(j); 

x(j)=xO+uO*(t(j)-tO)-0.5*g*(t(j)-tO) A 2; 

y(j)=sin(w*t(j)); 

v=w*cos(w*t(j)); 

uneg(i)=u0-g*(t(j)-t0); 

uplus(i)=2*v-uneg(i);u(j)=uplus(i); 

subplot(2,1,2);plot(thit(i),uplus(i),'o'); hold on 

u0=uplus(i);t0=t(j);x0=x(j); 

end 

m=size(x);n=m(1,2); 
sum=0;l=0; 
for k=2:n-1 
if u(k)<0 & u(k-1)>0 

amp=x(k-1)-u(k-1)*(x(k)-x(k-1))/(u(k)-u(k-1)); 
sum=sum+amp; 1=1+1; 
elseif u(k+1)<0 & u(k)==0 & u(k-1)>0; 

sum=sum+x(k);l=l+1; 
end 
end 

'The average amplitude of the ball is' 
sum/I 

subplot(2,1,1 );plot(t,x); xlabel('t');ylabel('ball position') 
subplot(2,1,2);plot(t,u); xlabel('t');ylabel('ball velocity') 


title(' The points o are the points immediately after the hits) 
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_ y" CO/i 1 

'Problem 3-57 Solution For The Frequency Spectrum' 
clear 

xO=input('Enter the value of the initial position of the ball >0 : ') 

uO=input(’Enter the value of the initial velocity of the ball : ') 

N=input('Enter the desired number of hits : ') 
w1=input('Enter the lower bound of the table frequency w : ') 

w2=input('Enter the upper bound of the table frequency w : ') 

iw=input('Enter the desired increment of the table frequency w : ') 
'Wait..., this might take a long time' 
x00=x0;u00=u0; 
nw=0; 

for w=w1:iw:w2 
nw=nw+1 ;wf(nw)=w; 
clear t x u y 
inc=(2*pi)/(300*w); 
g=9.81; 

x0=x00;u0=u00; 

t0=0;t(1)=0;x(1 )=x0;u(1 )=u0;y(1 )=0; 

j=i; 

for i=1:N 

while x(j)>y(j) I x(j)==y(j) 

j=j+i; 

t(j)=t(j-1)+inc; 

x(j)=x0+u0*(t(j)-t0)-0.5*g*(t(j)-t0) A 2; 

y(j)=sin(w*t(j)); 

u(j)=u0-g*(t(j)-t0); 

end 

t(j)=t(j-1 )+(y(j-1 )-x(j -1 ))*(t(j)-t(j-1 ))/((x(j)-y(j))-(x(j-1 )-y(j-1))); 
thit(i)=t(j); 

x(j)=x0+u0*(t(j)-t0)-0.5*g*(t(j)-t0) A 2; 

y(j)=sin(w*t(j)); 

v=w*cos(w*t(j)); 

uneg(i)=u0-g*(t(j)-t0); 

upius(i)=2*v-uneg(i);u(j)=uplus(i); 

u0=uplus(i);t0=t(j);x0=x(j); 

end 

m=size(x);n=m(1,2); 
sum=0;l=0; 
for k=2:n-1 
if u(k)<0 & u(k-1)>0 

amp=x(k-1 )-u(k-1 )*(x(k)-x(k-1 ))/(u(k)-u(k-1)); 
sum=sum+amp; 1=1+1; 
elseif u(k+1)<0 & u(k)==0 & u(k-1)>0; 
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sum=sum+x(k);l=l+1; 
end 
end 

avamp(nw)=sum/l; 

end 

plot(wf,avamp);xlabel('Tab!e frequency w');ylabel('Amplitude average'); 
title('Problem 3-57 Solution For The Frequency Spectrum') 
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From Bcj. -4,-4 — V Hie response dun ne^ fjrst 

Step function is 
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(U) = CO* ‘SjovaJ 1 itj 

' !q ' / 

i chicly hex*me Jnrfja.1 C-encfj-f/tns A** fye secondf?h*se , 
5/nce. <Qxcj fa?/0/) /S jtero d-nr/nj Me t/me A>t /J Ahe 
d-bore eyuatrons < 2 fc f/ftCfed /otb Me ‘feneraj solution 

For Me rrcc. (// braf/on 

?Ft) - Sftuu t -t- Ocft,) OrtcoJ: -t>t, 

Co)— ** ' 
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4-A 


H//'M zero t-mt/a./ Co nfAf/orji F<j. 4.2- / 

for fAe part/ cu /a r infejraj ,t4e e^uatto/j 


I 

*■(*) *J 


where f(f) '•* orh/'tarp exc/fetton P° KC -* . 

The S’ohti'/en f° r tAe Aomop&neou s epu&f/on 

■found Prom ~f> 

* 

y,(t) - Sprt,U> tM t -i~ c*5 oo^t 

and -til csmp/ete Solu 6/on As PAs Sum of these, 
ezju ait onf . <£ 

%,({)= S*rucd M t + * (o)C<va)J + J/(fj Scn.&jft-Jjff 

SP ^ 

- — 


From Pej. 4.2- / , t/>e response to o. unj / Step 

Function !$ f 

$(*> = y 

o 

Differentiate Aolh Sides «t. r. t. time 

ft*) = %(*/ 
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E<,. 4 .x-I * = fevtct-w 
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Fr.* Tr.L.4-7 f(t-f) = < (t-Xj 
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= to uct) — I (TO & /, y.Ct) - to£(.t) — zero 
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<-u 
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&J 
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* ct)** f ct)-xv^Lt-W *' 

... jS+/c M icF. is C(t-t,) 

For S *<. V"™"" 


* «; = *c* ; + 
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- 4 — J 1 — FT t 
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* 

(je* Sd. V-(jt) m fyS^ujrt ■h3C&u> m t f- 


,, *.ir 
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ArjS The response. to the ramp -Qnct/on of 

F, s 4.4-t ,s e /yen h y £fA*i -3 *-* 

*#- f[ /_ ^ 7 mt * 

p/M're.ntUt'O] && ej«*&*'* *»d &}uat,ry ^ +° 

•z.cro "the. peak time tp is obtained <*-5 

N j. / — Cos eJ«t/ 

Tin ojT- p — -r——--s- — ' 

$/nc£. c^tp muit he greeter -tha.n i ue. aXso ohtmn 

SLncJjtp = - > f^~C /-0»aV<i ) 

cw, - 

fx(,~ Oz *J„t,) 

Su.yH'huh^ these quantifies into Deft) fa/«** amfatude 
tS fouln <t 4S /-- ~r- 

&) - ' + ^7 1 t( " ^ 

^ r o'>r»* 

Lett in j T= be the period «t the 6Sc,Hator t the aJ>o*z. e^uetien is 

ptotted O-jT.nSt t'Ay ** '* F '<Z V4 ~‘8. 
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r*t 

J s^T&Ar 

* 

^Response 
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* * *■ 

jfl s*^, ft>_ t ^ 
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^ F » J**** *)r 

ip m ^ 4 - 7Ll>£ response must occur 
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II 

'U 

i 

r+- 

IV 


(ti \ _ / - xt«£ + *.rr«(S-tj\ 

(v “(,-i-^r 
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* 

<v 
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= ). F4- 


cj = i2C = 
* .ro 
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= /, * r * ^ F = 1.0 Kr.S. - lx.40* F„ 

.'. F« ‘ 10 
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^——— 2 4 —JZ 

% = |£- j A «« <&,(£-*.?) - ~ -,..j - « ^ j = o 

*' z ^ y^'(£ -oso)-**t x l£ i (^ ~’J - e *’^ i r' ~ ° 

\fjkich C ^y 


1"_ V//4 SufftfStJ Sm.hft Ar >f C*>f am.U A S+*. $> 

* =• ^-| <±M4*coj: — c*r>*±f s**><£ | s 4. —4>) 

Al#* rtsf*** *cc (*■'**-+) = f 

* 7 wkttt AsCO ” /[ ficS)vLuj d.if+{(&)<*»£Jjf 

0 0 


A -2: 


Ui £=+, *,« ± = S± = * 
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b 
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Caspt* - =■ '-ZS^l('f>+*±)t m * (f>-co m )t^O 
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"Therefore Ike, i/nte. t», (Lorrespond me 60 (% 2 ) <s 
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t - JL*2L 

071 p+^M, 


Su.bstifut.mq mto ~ksL we. hare 
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(M)^ -" -rh' r *<•&£ )-*M ] 
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^_part (b) for t > t| 

W)*ik fke -T ^<-co.fioir\ -4.^ — (£-) £ 



\/j\\Cr& t 
^Kafe to 


— A' f , 
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The r.rst fern cen fe reduced as &>/»«* 

sdJm&i.) = s~ *«*'(**> -o - -20Z ) 

\ / ~t~/\s J l + fir 

= S+H.z*ir 
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S^J 2 ,K ' 7r izk l )'\ 
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/.OSS 

3.0 

/./70 


88 



= sir - *s e (u(t) - t/t t J 



/4cc e/ jt 

a. = - i^J = 

$<V 

.5 «.£yzf 

J 0 JzG 

4.Z-S 



z 

... _ 15- f * 

" ^ / ff<X> 

c? 


o < i < i t 


= IS 

^ — S<^ < 

+£t,0-‘*>“~ t ) } 

0 <~t < t, 

d-7. 

d£ 

II 

fi? 

— OJ^Cmu^t + -jf Sj- — ° 


< 



"" t*nu> m t r ~ Cot / 

for 

\ 


“i,*, 

S4*C <X^„ = — 
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c '* eo ~ t r ~ '( J + M) i - 



/. Z<A* f j_ _ / _ - ] 

~ 1 w m £ # “*ihToZtf y if m,)'’ J 


4--Z5 


Jf~ t >'t / //e So/u./ton /S 

itrj-ij 


fr 

J L < 
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Z + Z = - — /<►© -XcS(t) 

, 12T s <7. £zr *. 

G*Jyoo$G £ki* — ^ # O& ^ 

Fir Me impose response. f it a* <x<*V/\s *JtJe j* - # treat *t separate/y. 

Ifs response /S - l£. ^S(f)S**.U)Jt-S)<lf = -lS+i*)0t 

e> 


Remain/ny So/ufto n Can !><■ CrOmpn.Fiij nu. mert Ca,//y -£rom D.£ 


CA.) 

(*>) 


■ - /tro — t w z 


z.-fit = 


7 v | A/V/i/ Con*. 

V-/ ~* 2 ; ■*■ z *>/ J a,-*/ = ° 


Siqrt *iM 21 j 

From E<f+J~f) 

-Fr* »*(*.) Z x 

Find Z ; front 


S. J (TO f **o ryy C**) 

= £ (o.osf tov = O./zso 

- /(TO - /ev (e.llsv) 

£j. ck) tie . 


H)hn\aHL cna G^nicl u$(- but IF tn o£>/■■ s. c/uc>' /'I 

n hx^ah. 


%This is the code for problan 4-26 

t0=0; 

tf=.05; 

x0=[0 0 ]'; 

[t,x]=ode23 ('prob426' ,t0, tf ,x0) ; 
plot (t,x) 


function xdot=prob426 (t,x) 
xdot=[100-100*x(2);x(l)] ; 


90 



(' This is exairple 4.7-2 on page 113 in the text’) 

clear; 

clg; 

J=0; 

T=0; 

X=0; 

Xdd=0; 

F=0; 

dt=0.05; 

T(1)=0; 

F(1)=0; 

X(1)=0; 

J(l)=l; 
for 1=2:25, 

J(I)=I; 

T(I) =T(1) +dt* (1-1) ; 
if I > 2 
if I < 5 

F(I)=500*dt*(1-1); 
else 

if I < 9 

F(I)=200-500*dt*(1-1) ; 
else 

F(I)=0; 

end 

end 

X(I)=Xdd (1-1)*dt~2-X (1-2)+2 *X(1-1); 

Xdd(I)=2*F(I)-16*pi~2*X(I); 
else 

F (I)=500*dt*(1-1); 

X(I) =dt /S 2*F(I) / (3+8*pi^2*dt^2) ; 

Xdd(I)=2*F(I) -16*pi^2*X(I) ; 
end 
ard 

A=[' J ' ' TIME ' ' DISPL ' 'ACCLRIN 

B=[J' T' X' Xdd' F'] 

plot (T,X) 

hold on 

plot (T,X, 'x') 

xlabel( 1 Time, seconds') 

ylabel ('Diplacement X, m') 

hold off 


'FORCE' 
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m a. 




;< ev/Jud/t* from <-0 


c*o 

cj; 


_3Toi,t . ,_ , / . /v on&t 

*ct) =x, e *~w(/7=p<** +*,)- -^3r 
* -s’**.* 

*ft; - J r )' r7 *** 


"Result *s .' 
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%Ct)=X t €. S***(h-s'e^t +<fi / ) = X& S*yf.(c^t-h^>) 

&r J* S'nO.ft 
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= Xfl^onj, = x t *>,S~4> = ° 

.'. <4> - tfO* a.net )( — M*. 

V 7 / 


7^£>7 7^2 ms. So/. H/i/A a-t>ovt i'lt'tii.J Conet. '■* 

p- 

At r^-t » sr /.o 

•*«« Pfek4-3Z 
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in the text') 


% ( 1 This is exanple 4.7-1 

clear; clg; 

J=0; 

T=0; 

X=0; 

dt=0.02; 

X (1)=0; 

Xd(l)=0; 
for 1=1:24; 

J(I)=I; 

T(I)=dt* (1-1) ; 
if I < 6 
F(I)=100; 
elseif I < 11; 

F (I) =-1000* (T(I) -0.10)+100; 


F(I)=0; 

end 

Xdd(I)=0.25*F(I)-500*X(I) ; 
if I > 1 

X(I+1) =2*X(I) -X(I-l) +dt /v 2*Xdd(I); 
else 

X (2) =X (1) 4<lt*Xd (1)+0.5 *dt^2 *Xdi (1) ; 

end 

axl 

J(I+1)=1+1; 

T(I+l)=dt*I; 

[' J ’ ' Time 1 ' Diplacanent ’ ] 

[J' T' X'] 

plot (T,X, 'O') 

hold cn 

tl=0:0.02:0.1; 

xl=0.05*(1-cos(22.36*tl)); 

plot(tl,xl) 

hold cn 

t2=0.1:0.02:0.2; 

x2=0.05*(1-cos(22.36*t2))-(0.5*(t2-0.1)-0.02236*sin(22.36*(t2-0.10))); 
plot(t2,x2) 
hold cn 

t3=0.2:0.02:0.48; 

x3=0.05*(1-cos(22.36*t3))-(0.5*(t3-0.1)-0.02236*sin (22.36*(t3-0.10)))+0.5*(t3-0 

2)-0.02236*sin(22.36*(t3-0.20)); 

plot(t3,x3) 

xlabeK't, seconds') 

ylabel ( 1 Diplacanent, m') 

hold off 
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%Ihis is the file f.m for problem 4-41 
function [force] = f(t). 
if t<.2 
force=500*t; 
elseif t<.4 
force=200-500*t; 
else 
force=0; 
end 
and 


problem 4-41 



1 




%T l his is the function for 4-42 
function [force] = f(t) 
if t<.2 
force=50*t; 
elseif t<-4 
force=10; 
elseif t<.8 
force=25*(t-.8); 
else 
force=0; 
ard 
ad 


problem 4-42a 




il-ij'X c.o/1 /■" 

%Ihis is t±ie function for 4-42b 
function [force] = f-(t) 
if t<l 

force=10*sin(pi*t); 

else 

force=0; 

aid 

and 


problem 4-42b 




4--4Q 

The stiffness rtx- cr 4'*z boom is represented Ay 
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%Hhis is the function for 4-44 
function [force] = f(t) 
if tc.l 
force=100; 
elseif t<.2 
force=1000*t-200; 
else 

force -0; 

and 

and 


problem 4-44 
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Displacement 
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clear 

global m k a step r 
‘Problem 4-46 Solution’ 
xO (1)=0;xO(2)=0;10=0; 

rreinput (' Enter the value of the mass (kg) : ') 

k= input (' Enter the value of the spring stiffness : ') 

a=input ( 1 Enter the value of the dartping coefficient a : 1 ) 

choice=input('Enter your choice of the forcing function function step=l and 

ranp=2 : ') 

if choice==l 

step=input (' Enter the value of the step function : ’) 
tf=input (' Enter the value of the final tine : 1 ) 

[t,x]=ode45('p449s',t0,tf,x0); 

plot (t,x(: ,1)) ;xlabel (’time') ;ylabel(’position’);title('Problem 4-4 :Step 
function response’) 
end 

if choice==2 

r=input (' Enter the slcpe of the ranp function : ') 
tf=input (' Enter the value of the final tine : ') 

[t,x]=ode45( 'p449r' ,t0,tf,x0) ; 

plot (t,x(: ,1)) ;xlabel('time') ;ylabel('position');title('Problem 4-4 :Rarrp 
function response') 
and 

function xdot=p449r(t,x) 
global m k a r 
if x(2)>0 I x(2)==0 
xdot=[x(2); 

(1/m) * (-a*x (2) ~2-k*x (1) +r*t) ] ; 

else 

xdot= [x(2) ; 

(1/m) * (a*x(2) ^2-k*x(l) +r*t) ] ; 
and 


function xdot=p449s(t,x) 
global m k a step 
if x(2)>0 | x(2)==0 
xdot=[x(2); 

(1/m) * (-a*x (2) /N 2-k*x(l)+step) ] ; 

else 

xdot=[x(2); 

(1/m) * (a*x(2)~2-k*x(l)+step) ] ; 
and 


109 



position 


Lj - £0/1 t— 


Problem 4-46:Ramp function response 
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position 
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clear 

global m k c mu step r 
'Problem 4-47 Solution' 
x0(1 )=0;x0(2)=0;t0=0; 

m=input('Enter the value of the mass (kg) : ') 
c=input('Enter the value of the damping coefficient (N.s/m) : ') 
k=input('Enter the value of the spring stiffness (N/m): ') 
mu=input('Enter the value of mu (N/m A 3) : ') 

choice=input('Enter your choice of the forcing function function step=1 and ramp=2 : ') 
if choice==1 

step=input('Enter the value of the step function : ') 
tf=input('Enter the value of the final time : ') 

[t,x]=ode45('p447s',t0,tf,x0); 

plot(t,x(:,1 ));xlabel('time');ylabel('position');title('Problem 4-47:Step function 
response') 
end 

if choice==2 

r=input('Enter the slope of the ramp function : ') 
tf=input('Enter the value of the final time : ') 

[t,x]=ode45('p447r',t0,tf,x0); 

plot(t,x(:,1 ));xlabel('time');ylabel('position');title('Problem 4-47:Ramp function 
response') 
end 


function xdot=p447r(t,x) 
global m k c mu r 
xdot=[x(2); 

(1/m)*(-c*x(2)-k*x(1 )+mu*x(1) A 3+r*t)]; 


function xdot=p447s(t,x) 
global m k c mu step 
xdot=[x(2); 

(1/m)*(-c*x(2)-k*x(1 )+mu*x(1) A 3+step)]; 
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clear 

global m k c mu step r 
'Problem 4-47 Solution' 
x0(1 )=0;x0(2)=0;t0=0; 

m=input('Enter the value of the mass (kg) : ') 
c=input('Enter the value of the damping coefficient (N.s/m) : ') 
k=input('Enter the value of the spring stiffness (N/m): ') 
mu=input('Enter the value of mu (N/m A 3) : ') 

choice=input('Enter your choice of the forcing function function step=1 and ramp=2 
if choice==1 

step=input('Enter the value of the step function : ') 
tf=input('Enter the value of the final time : ') 

[t,x]=ode45('p447s',tO,tf,xO); 

plot(t,x(:,1 ));xlabel('time');ylabel('position');title('Problem 4-47:Step function 
response') 
end 

if choice==2 

r=input('Enter the slope of the ramp function : ') 
tf=input('Enter the value of the final time : ') 

[t,x]=ode45('p447r',tO,tf,xO); 

plot(t,x(:,1 ));xlabel('time');ylabel('position');title('Problem 4-47:Ramp function 
response') 
end 


function xdot=p447r(t,x) 
global m k c mu r 
xdot=[x(2); 

(1/m)*(-c*x(2)-k*x(1 )+mu*x(1 ) A 3+r*t)]; 


function xdot=p447s(t,x) 
global m k c mu step 
xdot=[x(2); 

(1/m)*(-c*x(2)-k*x(1 )+mu*x(1 ) A 3+step)]; 
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clear 

global mu step r 
'Problem 4-48 Solution' 
x0(1 )=0;x0(2)=0;t0=0; 
mu=input('Enter the value mu : ') 

choice=input('Enter your choice of the forcing function function step=1 and ramp=2 : ') 
if choice==1 

step=input('Enter the value of the step function : ') 
tf=input('Enter the value of the final time : ') 

[t,x]=ode45('p448s',t0,tf,x0); 

plot(t,x(:, 1));xlabel('time');ylabel('position');title('Problem 4-48:Step function 
response') 
end 

if choice==2 

r=input('Enter the slope of the ramp function : ') 
tf=input('Enter the value of the final time : ') 

[t,x]=ode45('p448r',t0,tf,x0); 

plot(t,x(:,1));xlabel('time');ylabel('position');title('Problem 4-48:Ramp function 
response') 
end 


function xdot=p448s(t,x) 
global mu step 
xdot=[x(2); 

mu*x(2)*(1 -x(1 ) A 2)-x(1 )+step]; 


function xdot=p448r(t,x) 
global mu r 
xdot=[x(2); 

mu*x( 2 )*( 1 -x( 1 ) A 2 )-x( 1 )+r*t]; 
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£--38 
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Fron i E<f. 5*. J-S 7^e » O-l^rcJ 
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s-sz 


y \ —vw —([ — A >vv | —EED-AW--J' 


V_1 -4 A 

* - * 


-i— A4 + • t~ 


= * * 


H = «.?•?/ /f 

e *«*.** f£ 


Lit 4 — *”■—/, tUn y t =1 %. 35 s*A t y ~ 2.07 s«t. 

Lit At = 0.10 „ ** ,+x. ^ 

= *c £* + ***• - *;-! 

Vu> = & At ‘ + *?* ■ ^-/ 

lnifta.1 Contt. (T-/) ^ ^ 

*(') - $ ( ') — %(*) ~ ° ) */ ( o = 3 

Frany Fcj.^.^S) X (X) - j- (Q, XO)' [jJQ \ =O.LU> 

Fro no Fcj (4. 1 2j U n d T^S •£,o r 

W 9 -) = i ^ ?(*-) = ~ 4 #<»J 

%(*-)[j + £ &t v J = x a£*x /. ty(n)-0.00437. 

Computer pro$r«.v*\ a*cf plot o£ results -/ol(ow/. 
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3 - ooa f 

%This is the program for problem 5-52 

clear;clg; 

dt=0.20; 

f=100 

x(l)=0; 

y(l)=0; 

dy(l)=0.0; 
dx(l)=f/3; 
t (1)=0; 

x(2)=0.5*dx(l)*dt /s 2; 

y(2) = (x(2)*dt /s 2)/ (6+4*dt /v 2) ; 

dx(2) = (-2*x(2)+y(2)+f)/3; 

dy(2)=x(2)-4*y(2); 

t(2)=dt 

for i=3:150 

if i<95 

f=100; 

else 

f=0; 

ail 

t(i)=(i-l)*dt; 

x (i) =dx (i-1) *dt"2+2*x (i-1) -x (i-2) ; 
y (i)=dy(i-l) *dt"2+2*y (i-1) -y (i-2); 
dx(i) = (-2*x(i)+y(i)+f)/3; 
dy(i)=x(i)-4*y(i); 
ail 


problem 5-52 
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5-53 


(OV £ =-4W 

j = 4-x/a (*-%) ~6*/o -/ ) 


.'. £ - - 4o;t *• ^43 


h*+Aio* 
(ifl* * 


V - 40* - m2- + (> 02 , 2 = tos^vt ' x 

^5* IX St«rf (la *■*/><*?*. f/o n 

X(x) - j- *(*.)&, - A [-*0 xctjttoyctj] 

_/2 «2 

, ®t iftv = -r f f oxty-*‘ m i>(*)’+to ?(*j] 

t t ' J 

, 6 [/ * '■£ £} - yfc) & £ + to £$d) 

K*.) - ~ ^ -—--- £*(*.) 4^2 ^ r fo tSi(x) ] 

to to 

*(*-) = 2.4cro ‘Qfr) / (j£ + ?40 + Z4e-o ) ) 

9 7 ' 1 St«rj, 

yd) -(/+ 4z C)W)/ 4a -£* 
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^-5- "3 Co /»r 

%This is the mat lab code for problan 5-53 

clear; clg; 

dt=0.05; 

f=100; 

x(l)=0.0; 

y(l)=0.0; 

c^(l)=0.0; 

dx(l)=0.0; 

t(l)=0; 

z(l)=0; 

t (2)=dt; 

z (2) =10*sin (pi*dt) ; 

x(2) = (400*6*dt /N 4*z (2)) / (36+840*dt /s 2+2400*dt / '4) ; 

y(2) = (l+40/6*dt / '2)*x(2)/ (40/6*dt~2) ; 

dx(2)=-40*x(2)+40*y (2); 

dy(2)=40*x(2)-100*y (2)+60*z(2); 

for i=3:200 

t (i) = (i—1) *dt; 

if t(i)<4 

z(i)=10*sin(pi*t(i)); 

else 

z(i)=0; 

end 

x(i)=dx(i-l) *dt /s 2+2*x(i-l)-x(i-2) ; 
y(i)=dy(i-l) *dt' x 2+2*y(i-l)-y(i-2) ; 
dx(i)=-40*x(i)+40*y(i); 
dy(i)=40*x(i)-100*y(i)+60*z(i) ; 
end 
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wJ & } - -rniJ$(/!&! - COf - dt) 
-rnds,'/i9i t k&si'/fCQf-fy) 


clear 

global m I a k g 

'Problem 5-63 Solution' 

m=input('Enter the value of the mass (kg) : ') 

k=input('Enter the value of the spring stiffness (N/m): ') 

l=input('Enter the value of the pendulum length (m) : ') 

a=input('Enter the value of the distance a (m): ') 

xO=input('Enter the four dimensional initial condition [xl xld x2 x2d] : ') 
tf=input('Enter the value of the final time (s) : ') 
t0=0;g=9.81; 

[t,x]=ode45('p563d',t0,tf,x0); 

subplot(2,1,1 );plot(t,x(:,1 ));xlabel('time , );ylabel('theta_1 ');tit!e('Problem 5-63') 
subplot(2,1,2);plot(t,x(:,3));xlabel('time');ylabel('theta_2'); 


function xdot=p563d(t,x) 
global m I a k g 
xdot=[x(2); 

-(g/l)*sin(x(1 ))-(k/m)*(a/l) A 2*sin(x(1)-x(3)); 

x(4); 

-(g/l)*sin(x(3))+(k/m)*(a/l) A 2*sin(x(1)-x(3))]; 


149 




I~ ( ( X, J Xt) — (x t Xf J — f[' C°y X±- X£) J 

I- y= X z -x, .Therefore, 
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m fe 
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w P ,/•/ > ./« ' a.Ccel(rat lesn 
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Ip the 
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J/'rt ebrort 
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Wt 



"*1= 7"P'1* S a t fs the. tapje/rtto/ dccele^at/oii) 
m C 16 - xsin &) = 

" 16 - X sin 9 f pin 9—0 

The -1 \Uo.tfons of notion 
j< - 19 zipd-/^cos9 ffe)A = 0 

16 - x sine + pin 9 - o 1 


To fi»d tht natural freft/nci'es cue a-sut^re 

X = X linu/t and &=9 sipu/t " ' 00 

It r> clear- {hat 0=0 y (M- fk/m 1 
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%This is the function file 
function [force]=f(t) 
if t < 0.1 

force = 10* t; 
elseif t < 0.3 
force=l-(1.7/.2)*(t-.l); 
elseif t<0.5 

force=-.7+(1.2/.2)*(t—.3) 
elseif t<0.7 
force=.5-4*(t-.5); 
elseif t<0.9 
force=-.3+(5/2)*(t-.7); 
elseif t<l 

force=.2-2*(t-.9); 
else 
force=0; 
and 
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clear 

'Problem 7-26 Solution' 
global m k 

m=input('Enter the value of the mass m (kg) : ') 
k=input('Enter the value of the spring constant k (N/m) : ') 

xO=input('Enter the value of the four dimensional initial condition [x xd theta thetad] . 
') 

tf=input('Enter the value of the final time : ) 
t0=0; 

[t,x]=ode45('p726d',t0,tf,x0); 

subplot(2,1,1 );plot(t,x(:,1 ));xlabel('time l );ylabel('x');title('Problem 7-26 ) 
subplot(2,1,2) ;plot(t,x(:,3));xlabel ('time') ;y label ('theta'); 


function xdot=p726d(t,x) 
global m k 
xdot=[x(2); 

-4*(k/m)*x(1); 

x (4); 

0 ]; 
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clear 

'Problem 7-27 Solution' 
global m k 

m=input('Enter the value of the mass m (kg) : ') 
k=input('Enter the value of the spring constant k (N/m) : ') 

xO=input('Enter the value of the four dimensional initial condition [x2 x2d x3 x3d] : ') 

tf=input('Enter the value of the final time : ') 

t0=0; 

[t,x]=ode45('p727d',t0,tf,x0); 
z=(1/3)*(x(:,1 )+x(:,3)); 

subplot(3,1,1 );plot(t,z);xlabel('time');ylabel('x1') ;title(' Problem 7-27') 
subplot(3,1,2);plot(t,x(:,1 ));xlabel('time');ylabel('x2'); 
subplot(3,1,3);plot(t,x(:,3));xlabel('time');ylabel('x3'); 


function xdot=p727d(t,x) 
global m k 
xdot=[x( 2 ); 

-(k/(3*m))*(2*x(1)-x(3)); 

■ x(4); 

-(k/(3*m))*(2*x(3)-x(1))]; 
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>sl v i= l z u u ; u x u ; u u xj 

S'-3. CO'jyL 
M = 

2 0 0 • 

0 10 

0 0 1 

»K=[ 3 -1 0 ; -1 2 -1; 0-11] 

K = 

3-10 
-1 2 -1 

0-11 

»t=inv(M) *K 


t - 


1.5000 
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0 


-0.5000 

2.0000 
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0 

- 1.0000 

1.0000 


»[V, D] =eig (t) 

V " <t>, <h <h 

0.7569 0.3031 0.2333 

0.2189 -0.8422 0.5808 

-0.6158 0.4458 0.7799 


D = 


th ese- can 



/, S'Sr9^~ i 

/ . 0 0 a 


1.3554 0 

0 2.8892 

0 0 


»croega=sqrt (D) ^ 



0 

0 

0.2554 


crnega = 

1.1642 0 0 

0 1.6998 0 

0 0 0.5053 
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Lzl 


CM or* cterts //<• equation in "terms of stiffness 


\ - f 0 

-J (.*->;) 'I 
c -> (i-M 


- o 




o\ C J 

A - 4> 


mode Z S utsstttintd *X — J • 3 


•HU -/ o 

-/ .6444 

o -/ -.3i34 


'«,] f° ] 

u l 5 id I 

u;i Uj 




■p /-* 'ZS^Z and ccdd to l ’"'rod 


l - 2.+87? 0 

0 -Z.VI3Z 

L o -/ 

nt 



- 0 


From 3~ oltdif) 

» l S - ro»s % f -.35S+(- 9.4*7*) zO -I.MSj 

. d> J' AZ ^l 

l /. j 

. (• .(>800 | 

For <^b do Sdmt, to obtain V 

j J ] ,.m> I 

Frti $-1 comp«ter results . 


i -* 


f(3 - Z*) “ 1 0 

-i (*.- >; -/ 

O _/ (>-» 


“ ° bntko^ o-P ^factors // 

is not necessary to eJter +b<- 


, A x -v 

~M (•- 


C a= “ 


>) 


-l -l 

o (.'-*) 


fi 'rsi ro'd 

ut \ t - .l&'+ C » ~ 

C.*- *7447 
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5? — ± C*n±. 


C 'S ~ 


-/ (*->) 

O -/ 


= 1,0 


•LI? I 
4> = t -7447 


/. 00 


Fpr 4 S«£ A x = etc . into Stye CoP'actcrs 

Can a./so Sta't WitJi P'/cta.l/j/fy aj-eftcn 


t-5 


*P met/on m terms oP CoPhcfcrs o"P a. ce/umn >S 
Dot prcyorPtene-I ^ e-f^en \/e.ctors . 

t-C _ 


h$f' ,Y E ~i °, ] 

Gam$ kl\m . ^X, = 


)-tci -i 0 
-I t-rfl -/ 
o -) .scjl 

AjLC 

■■■ 4. *| .m 

1.000 


~\ z * !'?S£ 


(*-*> -' 0 
- -/ (a-X) -1 

O -/ ('-*> 


^ o 


/\ U>J»* 

A = 



* = *So2 

z 


= 0 


% - i^ — - . 44 i" 

' / .*01 

whtn +k«r£ AOt XtnpS in mtf.fYn> 
Ga>tS '»S net necessary 

r- o. 3 > 3 p re b )c »n 


’ .44^ 

-) 0 1 

fx. 1 ) 

^ = 1.0 

II 


-1 

.445 -/ 

l*l( = 0 



o 

-i -.ssr 

M 

Xf - -7.14] 

( 

1.000 


X 3 - Z.1A1 


- I.X47 -I 0 (*' j 

~i i+i y* > z0 

O -/ -2.I47JJ/.0 


^-/.o 
*, : -2M7 

X r /.&?! 


4- 
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•g'-6 Cc/i'K 

» D =[ 2 -1 0 ; -1 2 -1 ;0 -1 1 ] 
D = 


2 -1 

0 


-1 2 

-1 


o 

1 

1 


»[V,D]=eig(D) 



V = 





<£3 

0.7370 - 

-0.3280 

0.5910 

0.3280 - 

-0.5910 

-0.7370 

-0.5910 - 

-0.7370 

0.3280 

D = 



1.5550 

0 

0 

0 

0.1981 

0 

0 

0 

3.2470 

»sqrt (D) ^ u3 

f IT 


ans = 



1.2470 

0 

0 

0 

0.4450 

0 

0 

0 

1.8019 


These. 


K) 0> ) ted 

,VVJ-OI 




/, 0 0 o 

/-/.VO'S-\ 
= -3.a^ I 

1 / .t O O J 


» 
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{?r iX 


From Example £>./-/, +** f/e»il;h‘iy matn 


/S 


l 

k 


/ / ' 
/ 2 2 
/ 2 3 




ha eauo.1t in o/ molten 'S 


[; i f] -&~)i 

J 2 5 J 

Charoxter/sft O' *(jU*hon 


- -*» 

Wj 


=. O 


/.ef > = 

r r> 


=. o 


(»->> ' ' 

/ U->) * 

f a C3 - a) 

> -S.i>SbS 

- 4.4$Z>.r / 7 

/ -J.DSOS Z 
, a —2,osoiT 

Odd ro«!S I andi to 6r"> '**' reyJ 1 


A, - OS-OS 
*,.= 0.445 / 
^ - O.JO^O 


/*,»«/ ity 4. OSOS 


-I .246<j -ZACI 
/ -J.A&T 2 
/ 2 -a. Graft 


- 4.CS0S ! 1 

1 Yrttf cJi ro*J$/**} *£• form r^u/5 

Q -Z.WL *- 24c 1 

/ 2. - J ^ rar 

1 

| 0 Z.iAll -I-Z034. 


/'ejt/ /i* '* 

- 4.0SOS 

o 

0 


/ / 

-2-ML 2.24i<j 

£ 2 . 44 ? -I.&SC 


C*Mt,n«e te reduce term 5> 
t5> 


% = *•***! = .&>/* 

1 Jt.?054 

^ . i-8t>r*L - .4447 

*i " 4.<?^r 


~4.tS6S l 1 

0 -XX01L i.14^ 

O © 


ft, 

J'**. 


D‘. 


444 7 ' 


O 

See^Bfcb £-4 


0 ; 


Uo 


4,=j- f<, ' + > £, 4 , * A 

J . £><?0 I 2 ^ 


/ /-/.247 1 f/-S/? 

\ = 0^}| «»<< V •***«> !')!« 
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% ~ $ £4 4*4 '*»• per 4sssy/>e4 bp jnsArttcfcr- 

■ g — 3 V/J# ^ a** « *rt^ 7^ £ijU&/r*t) of 

* 

*7161 ton IS 


3 -J - i 0 

1 



’ 0 * 

-/ Z - / O 

\ 

tl *,\ 

• 4 

0 

-/ -/ J -/ 

! 

] */ 1 

* — 

1°, 

0 O - ! ! 


' L'X+I 


i 0 J 


A to <VH 

A ' "X 


CM4racfertff/c eyuoiton 


(3-*) -/ -/ 0 

-/ (1-A) -/ 0 

-y -/ c3->; 

o o -i l'-x) 


i cj- v . „ 



(*-x) 0 


-/ -1 0 j 

-/ (Jt“X) 0 


=p-v 

-1 (3-sV-) 

+ 

-/ - 

-/ -/ -j 

= 0 


0 -/ 0-K 


0 -. H 

0 0 C>-X) 



-i -J 
6 ('-X) 


(3->) - J 

~ ) o-M 


][t5-*>(‘->H3 -[C7-iA)(*-x)j = o 
•’• "X A -5 “X s ■+ J.S* >* -*/ h + 3 =■ O 


-1 —J 

o 0~» 
+ (*-X) 


-/ -I 

o 0->> 


— i -/ 

O ('-X)| 


= o 


The. (Zioentfa/ucs Cat) wittier f J e. ^e/ernj/yj ee ^ Prom, 
roofs e>P V-<A£- (TKaraoheri^hc. e^u.<x^-r on or cAirechlp' 

U5;a, Re COtnnra^d Gift- . 


»c=[l -9 25 -21 3] 
c = 


1 ~9 25 -21 

»roots(c) 

ans = 

4.1796 

3.4882 

1.1538 

0.1783 
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CoM. 


3/5 


»I^[ 3-1-10 ; -12-10 ; -1 -1 3 -1 ; 0 0 -1 1] 
D = 

3-1-1 0 

-1 2-1 0 
-1 -1 3 -1 

0 0-11 

»[V,D]=eig(D) 


0.6371 -0.5800 0.3679 -0.3497 

-0.6679 -0.0883 0.5672 -0.4737 

0.3569 0.7725 0.1120 -0.5131 

-0.1434 -0.2430 -0.7283 -0.6245 


3.4882 000 

0 4.1796 0 0 

0 0 1.1538 0 

000 0.1783 
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fi - to 

Opt.Wl (S**- C*mp+er r ro^c*f*$ ftr?rtbi-^) 


% - if A DOF Aut© 



«\>£/ ' - * c */ - */[*/" 

= -M, -*,[*•- <*♦«»•>] 



ct,+^x) (^v^A) ~£) ~ 

bxtrt'Q (*A + V‘; *A 'M* 

- 4, *,£/ (*,+*.; o 

- Ax -V 1 - o (*,+**; 



jjzJrZ 4 DOF Auto 


Too | 


SC>t>0 

35"0^ 

1-2.40 0 

j 

-1UOO 

/60dl 

4,- 

3SD0 

127 250 
_/ _ 

^/O,soc 

- 14100 

pw _ 

Km^ m 

-2000 

! 0 , zoo 

| \4djtoO 

1 

0 

| 4/J^ 


-JJ40 

~/43 oo 

l 0 

4/000 


. 

— 
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If-fZ Co/if. 


»M=[ 100 0 0 0 ; 0 1600 0 0 ; 0 0 4.969 0 ; 0 0 0 4.969] 


M = 


1.0e+03 * 


0.1000 

0 

0 

0 


0 

1.6000 

0 

0 


0 

0 

0.0050 

0 


0 

0 

0 

0.0050 


»K= [ 5000 3500 -2400 -2600; 3500 127250 10800 -14300; 
-2400 10800 40800 0 ; -2600 -14300 0 41000] 


K = 


5000 

3500 

-2400 

-2600 

3500 

127250 

10800 

-14300 

-2400 

10800 

40800 

0 

-2600 

-14300 

0 

41000 


»D=inv (M)*K 


1.0e+03 * 

0.0500 0.0350 

0.0022 0.0795 

-0.4830 2.1735 

-0.5232 -2.8778 


-0.0240 -0.0260 

0.0067 -0.0089 

8.2109 0 

0 8.2512 


»[v,d]=eig(D) 


v - 


0.9940 

-0.0671 

0.0767 

0.0398 

f 0 4-tvW d/i 


-0.0031 

- 0.0011 

- 0.0202 

0.9998 


0.0030 

-0.0008 

-0.9998 

-0.0203 


10 

v i 0 A 


-0.6996 

-0.6473 

0.1314 

-0.2727 
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Coa+. 


1.0e+03 * 

ro^' 6A 

0.0447 

0 

0 

0 


UJ I , 

V i b r Q *i 0/\ 



8.2560 0 

0 8.2141 

0 0 


4^-045 /a^ 
0 
0 
0 

0.0768 


£ = 13 

I o itrv l 


[*] = — 

4* El 


/ 

-/.S' 




>>M=[500 0 ; 0 100] 
M = 

500 0 

0 100 


»K=[ 6 1.5 ; 1.5 1] 


K- = 


6.0000 1.5000 

1.5000 1.0000 

»d=inv(K)*M 

d = 

133.3333 -40.0000 

-200.0000 160.0000 


214 



is-—/3 Cont' 


12 & ez 


»[v,d] =eig (d) 


v = 


-0.4605 0.3597 

-0.8876 -0.9331 


d = 


56.2356 0 

0 237.0977 



*- l± 


©f mo'tien *S 



- 

< 

»• 

4*v 

2 


k 


1 


1*1 



»M=[ 3 0 0 ; 0 2 0 ; 001] 
M = 


3 

0 

0 

0 

2 

0 

0 

0 

1 

»K= [ 6 
— 

-2 0; 

-2 3 -1;0 

K = 



6 

-2 

0 

-2 

3 

-1 

0 

-1 

1 
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o O 



%~l U i COnf 

»d=inv(K)*M 


d = 


0.7500 

0.7500 

0.7500 

»[v,d]=eig(d) 


v = 


0.5000 

1.5000 

1.5000 


y 


0.2160 

0.5376 

0.8150 


-0.6256 

0.6743 

-0.3924 


0.2500 
0.7500 
1.7500 


0.4337 

0.3636 

-0.8245 


) 


d = 


2.9382 

0 

0 


0 

0.3678 

0 


0 

0 

0.6939 


u>, - . 

-hk>H9 

u;, - l.^oi ffy, 




±z IS 


M- 


TT* rn&lrix , s EUontf/e. 

Grtc. CryJL 


-Jl 

3El 


I* * 

Z1 14- + 

(l, \ 

14- 3 *.$" 

U 

4 X-S 1 

ll 

U -* 

K % j) 


J n V*.rst of &3o\re. nta.tr!} /$ 


3El 


‘ 27 

(4 

4 * 

-/ 


~ • *3* S' 

-1*3/ * 

14 

4 

s 

z.S 

rj 


_ 3*1 

“ 1 * 


3.12S -3.E3f 

L Jf4 
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CtA-h 

5 * 4 = [ 1 0 0 ; 0 1 0 ; 0 0 1 ] 


M = 

10 0 
0 10 

0 0 1 

»K= [ .5385 -1.231 .9231 ; -1.231 3.385 -3.539; .9231 -3.539 6.154] 


K = 


0.5385 -1.2310 

-1.2310 3.3850 

0.9231 -3.5390 

»d=inv (K) *M 


27.1245 14.0695 
14.0695 8.0387 
4.0223 2.5124 


0.9231 

-3.5390 

6.1540 




3_ex. 

CO z -JL s 


4.0223 

2.5124 

1.0040 


»[v,d]=eig(d) 


f, jx- 

0.4526 0.1737 0.8746 
-0.6822 -0.5642 0.4651 
-0.5743 0.8071 0.1369 



» 


0.8180 0 0 

0 0.1133 0 

0 0 35.2359 



OOV 
'2‘ <P 
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1 ? — 20 





K “ [-/ «] 

-/ _ T"» AA- 12 .| ./4-<?ar7*l 

L 0 


_T 

4472 0 “I{3 0 " 

0 M Q = 

. 1 ^/ .-*r+ 0 r 


[2.1*4/ -. 4411 ! £M 

Q s I O App CS 

fa T MQ''->T]Y«o 5-Ji, 


'.4472 

. i4<?r 


"".Coo 

.2.06 ' 

0 

•7454J 


• XOO 

r* 

•4222 



<i+*r\das4 Porn* £ yA — XI =. O V = e>r y : Q Y 


2 - 2 / 

k 


S«eA?p. C5 - 



-/ 

2 


- / 


-/ 


b' ^ 


''7? 2 -S774.-SJ-J+ 
° ‘‘*9/0 -I. 05 & 

0 0 -7^ 


cf 


I 



•S77$ 15?Z •llAC 

0 -7W6 l-OSiS 

0 0 *'WO 


»k=[ 3 -1 -1 7 -1 2 -1 ; -1 -1 2] 
k = 

3 -1 -1 

-1 2 -1 

-1 -1 2 

>>q=<±ol (k) 

q = 

1.7321 
0 
0 


-0.5774 -0.5774 

1.2910 -1.0328 

0 0.7746 


218 



g- Ifc 



»d=inv(k) 


3 2 1 
2 2 1 
111 


»[v,d]=eig(d) 


^ ^ A 

0.3280 -0.5910 0.7370 

-0.7370 0.3280 0.5910 

0.5910 0.7370 0.3280 


d = 




0.3080 

0 

0 


0 

0.6431 

0 


0 

0 

5.0489 



u->, - 


-/.avT'y 






'o 



g-17 






/ 



"a -/ <5 o" 

(*•' 

m 

/ 

/ 

7 * 

: + 1 

-» 1-/0 
o -/ i 

•»* 


/ 



-0 0-/2 




219 



~ l ~h C'OA^- 


»K=[ 2-100 

; -1 2 

-1-0 ; 0 -1 

2 -1; 0 

0 -1 2] 

K = 





T ^ 






2 

-1 

0 

0 



-1 

2 

-1 

0 



0 

-1 

2 

-1 



0 

0 

-1 

2 



>>d=inv(K) 






d = 






0.8000 


0.6000 

0.4000 

0.2000 


0.6000 


1.2000 

0.8000 

0.4000 


0.4000 


'0.8000 

1.2000 

0.6000 


0.2000 


0.4000 

0.6000 

0.8000 


»[v,d]=eig(d) 





v = 






0.3717 


0.6015 

0.6015 

-0.3717 


0.6015 


0.3717 

-0.3717 

0.6015 


0.6015 


-0.3717 

-0.3717 

-0.6015 


0.3717 


-0.6015 

0.6015 

0.3717 


d = 






2.6180 

0 


0 

0.7236 

0 

0 

0 

0 

-•-= M rr 

0 

0 


0 

0 

0.3820 

0 

0 

0.2764 

W <3 - 

- JTO ^ 
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't-'r 

»x=[2 -1 ; -1 4] 

X - 


2 -1 

-1 4 

»U=chol (x) 

U = 


1.4142 

0 

*U'*U 


-0.7071 

1.8708 

Check Cal do. loM 6^ 


ans = 

2.0000 - 1.0000 

-1.0000 4.0000 

r-iy 

»x=[3 -1 ; -1 1] 


X = 


3 -1 
-1 1 


»U=chol (x) 

U = 

1.7321 -0.5774 

0 0.8165 
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Co/it- 


»inv(g) 
ans = 


0.5774 0.2582 0.7746 

0 0.7746 1.0328 

0 0 1.2910 




»m=[3 0 ; 0 2] 


3 0 

0 2 

»k=[3 -1 ; -1 2] 
k = 

3 -1 

-1 2 

»q=sqrt (m) 

q = 

1.7321 0 

0 1.4142 
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£ Cnt 

»a=inv(g) *k*inv(q) 
a = 


1.0000 -0.4082 
-0.4082 1.0000 


»[v,d]=eig(a) 
v = 


^J2L 

K 


-0.70^1 -0.7071 

0.7071 -0.7071 




d = 


1.4082 

0 


0 

0.5918 


co, = 


S^-aj 


-mW 

—-OS - 0 
^ O /■ 


x*n= 


[4 0 0 


"(f 2 0"; 0 0 11 



5€C ^ycx/nplt 


m = 


4 0 

0 2 

0 o 


0 

0 

1 


»k=[4 -1 0 ; -1 2 -1 ; 0 -1 1] 
k = 

4-10 
-1 2 -1 
0-11 


D-e^c O'*!/? <5.5c 

K 
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'- - 

»a=chol(k) 


a = 

2.0000 - 0.5000 0 

0 1.3229 - 0.7559 

0 0 0.6547 


»ain=inv(a) 


ain = 

0.5000 0.1890 0.2182 

0 0.7559 0.8729 

0 0 1.5275 


»atild=ain' *m*ain 
atild = 

1.0000 0.3780 0.4364 

0.3780 1.2857 1.4846 

0.4364 1.4846 4.0476 

»[v, d] =eig (atild) 


v = 


0.4232 - 0.8944 

- 0.8516 - 0.3381 

0.3094 0.2928 


0.5585 0 

0 1.0000 

0 0 


0.1447 

0.4006 

0.9047 


A = 

0 * 

0 

4.7749 





OU, z /.XiSr/ ^ 

- /. 0 d 
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§~- £3 


>>x=X2m*v 


x = 


0.1182 -0.4472 0.3455 
-0.3737 -0.0000 1.0926 
0.4726 0.4472 1*3820 

»m=[4 1 0; 1 4 1; 012] 


4 10 
14 1 
0 12 

»k=[2 -1 0 ; -1 2 -1 ; 0 -1 1] 


x- o y 


/'e 

\_CC r rnQ~' 



Tke_ xyio^4f'i / \ proh/t'^x 
13 

”7T 

A ^ Q-y 


k = 


2-10 
-1 2 -1 
0-11 


>^=chol (m) 



2.0000 0.5000 0 

0 1.9365 0.5164 

0 0 1.3166 


225 



tf-M Co/it 


»qia=inv(<^) 
qin = 


0.5000 -0.1291 0.0506 

0 0.5164 -0.2025 

0 0 0.7596 


»a=qin' *k*qin 


a = 


a*>ic 


0.5000 

-0.3873 

0.1519 


-0.3873 

0.7000 

- 0.6668 


0.1519 

- 0.6668 

0.9923 


»[v,d]=eig(a) 


v = 


0.4887 0.8165 0.3074 
0.7480 -0.2108 -0.6293 
0.4490 -0.5375 0.7138 


d = 


0.0467 0 

0 0.5000 

0 0 

»x=qin*v 
x = 

0.1705 0.4082 

0.2953 -0.0000 

0.3410 -0.4082 


0 

0 

1.6456 


x= 

0.2711 

-0.4695 

0.5422 


» 
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[q~ r m a' - 7%] - £> 


- J<- 

^ " Uf-r*) 


X = gT7 


»k=[2 -1 0 ; -1 2 -1 ; 0 -1 1 ] 


2-10 
-1 2 -1 
0-11 

»m=[4 1 0 ; 1 4 1 ; 0 1 2] 


4 10 

14 1 

0 12 

»q=chol (k) 


1.4142 -0.7071 0 

0 1.2247 -0.8165 

0 0 0.5774 

»qin=inv (g) 

gin = 

0.7071 0.4082 0.5774 

0 0.8165 1.1547 

0 0 1.7321 

»a=gin 1 *m*gin 


2.0000 1.7321 2.4495 
1.7321 4.0000 7.0711 
2.4495 7.0711 18.0000 
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»[v,d]=eig(a) 


v = 


-0.8165 

-0.4714 

0.3333 


-0.1494 

-0.3851 

-0.9107 


0.5577 

-0.7934 

0.2440 


d = 


2.0000 0 

0 21.3923 

0 0 


»x=qin*v 



x = 


X - 



-0.5774 

0.0000 

0.5774 


-0.7887 

-1.3660 

-1.5774 


0.2113 

-0.3660 

0.4226 


%-Z(o 



c 


-tyx 


<H* f C 




*,x, 


Oft. 


n 


***** - - 



(AfK+h) " K " 
- \ (*S K) 


{* 



0 

o 
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C«/*r 


»m= [10; 01] 
m - 

1 0 

0 1 

»k=[ 3 -1 ; -1 2 ] 
k = 

3 -1 

-1 2 

»[v,d]=eig(k) 

v = 

-0.8507 -0.5257 

0.5257 -0.8507 


3.6180 0 ^ ' C ^? 1 

0 1.3820 K 
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; 

/ 

t , 



3 

- t 

-1 

0 


K = * 

-/ 

2 

0 

0 

3 

-1 

L 



O 

— 1 

-* 

z_ 


=. - 1c, + ^ (*,- */,) +■ ^ faj ~ */) 
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- 0.6018 

0 


ECU»eig(a) 


ans = 


0.3978 

0.0360 

1.5548 

0.5000 


- 0.0666 

0.4571 

0.3078 

0 


- 0.6018 

0.3078 

1.1350 

0 


0 

0 

0 

0.5000 



■ ISO -4-3Z i -.iso -.433" 

^ _-l_SO 1 — •433 -.ISO 

J 2 T - 2 ^ o ~ .433 


tleme^t a.-3 
<-£ 


=?3ov‘ 


C~ .5 
C A - - 25 o 


- Sh 

S*"- .7.T0 


cs = -.^23 




b. le e nf 3_£• 


276 



F ro*y\ E>a*rsp It (0,S“l <X*uA *f^£ £ rxcts l cx^cl 3 

are f')<ec( / +Ue. sY,f( ness ma.tr.> 4.Vf/accme«t Victory 

arc : 


r r i 


r — 

r — \ 

r t 


*4 O ! -12 3l 

f ^ 

M, 


o n 1 i -je -s* 1 

J«M 


r ej 

-12. -3£ | 12. “i£ 

A 

0 


) 

3i •5’/*’ j “ii £ l 


0 



M * = ? f[ 3£ ^ =*f[ 3 ^i) + - s£ ffi) 


10-15 


"Rcftr to 


I W 1 


= e - o 


_< JV £ f 3 12 0 


?40 


2 * 


SET 


24 ° 

o 2 e ; 



^ = e£r)(-4) = 

(24 -3IZ\) 


0 


oj, - li. 14 


0 

Il*( I - X) 


/frt £ £P 

6720 El 

1 

- O 


s 

6 T t 

e . 


^ ™ ^ 

' ~ 312. ~ 672.0 E] 


^* = 1 = 




Un £4 


cO t - Si -'!* 

tx«.ci vctliKj ( 


#■ 

u >, ^ f ~ 

oJ v = 6 /. 67 T 


-> t.C. °f D K^K 
-> SS^c H . 


6 720 ej 






10 - 17 



AA - **i 

M - SE 


K • 


~ /si. f-li 1 6*. —/;Si 

_ z _ zi . 

5'4 /// i ~* f % 

-/si -$ f\~ x * L A t l 


fZ 

U _ 

-/% 

a 


- 

fir, 


h 

< 

1*4 

J 

L *J 


a 

i ->1 

u ' 

4l % 

1 -U 

zi z 

-U , 

f/2 

-Ct 

zl* 

-6l 

4t z J 


Far fie f/rrxfe/-free. Seam 4^=0 u/htcj, e/,nn ra/tS fAe first 
Ca/umr , &*<£ fir the e.nfe*i proA/em, Fa /J J1 ° /y1 ^rest’ 
the firs/’ ra*J car aJso Se JeJetecL. Also u/Aer tAe cAerrsf/c. 

4qu action j -cfiW + K i = O /S SaSreet, tke 1 ,n tAe 
*r><Ltr/C& Factor out, dnd "*■ Can fit * * AO and 
te*/nt& the <ej as.‘ 




■4ZO 


A 

fZ 

-3 


/S -3 
/S-L -Z2, 
-ZZ A 



‘4 

—£> 

A 


-<o Z 
!Z -C 
-6 +. 


o 


// si,u!d kt acted. /had the stdVcess maf‘X Mas a. 

Setcrauaaai. «h ,cX /i zedc sad rfiere/ise car, net Ae 

,a verded cr dec.*, f,sad. It '■> tier,-fire accessary t, 

deem f>CSC tAi <W *,etns . 


279 



1 0 -ircOA,r 

4 13 -3 ; 13 156 -22; -3 -22 4] 
m = 


4 13 

-3 

13 156 

-22 

-3 -22 

4 

»k=[ 4 -6 2; 

-6 12 

k = 

4 -6 

2 

-6 12 

-6 

2 -6 

»q=chol (m) 

4 


q = 

2.0000 6.5000 -1.5000 

0 10.6654 -1.1486 

0 0 0.6563 

»Qin=iiiv' (q) 

gin = 

0.5000 -0.3047 0.6094 

0 0.0938 0.1641 

0 0 1.5236 

»a=qin 1 *k*qin 

a = 

1.0000 -0.8907 2.2503 

-0.8907 0.8198 -2.3868 

2.2503 -2.3868 10.6088 
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/O-'l?- <ZCA + . 


»[v,d]=eig(a) 
v = 

0.5916 0.7760 0.2186 
0.8043 -0.5495 -0.2262 
0.0555 -0.3096 0.9492 
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/O 


» 

What is the length of the bar? 1 
len = 


1 


How many elanents would you like? 2 
n = 


2 


What is the mass of the beam ? 1 
mass = 

1 


What is the elastic modulus of the beam ? 1 
E - 


1 


What is the mcment of inertia of the beam?l 
I = 

1 


v = 


0.1311 

0.4801 

-0.6584 

0.5245 

0.2092 

-0.0266 

-0.3339 

0.3855 

0.5244 

0.6813 

0.3884 

0.6812 

0.4446 

-0.2905 

0.3211 

-0.2459 

-0.1322 

-0.4223 

-0.5961 

0.6233 

0.8779 

-0.4202 

-0.2049 

-0.1008 

0.0219 
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/O-ttr C,or\r 

lambda = 

1.0e+04 * 

6.2116 0000 
0 1.8376 000 

0 0 0.3116 0 0 

000 0.0241 0 

0000 0.0000 


cmega = 

249.2306 0000 

0 135.5595 000 

0 0 55.8228 0 0 

000 15.5142 0 

0000 0.0000 

X = 


22.2515 

29.2961 

-13.3294 

5.4582 

1.7321 

0.4353 

-1.1285 

0.4848 

1.1871 

0.8660 

28.3733 

10.7305 

11.8711 

-2.7048 

1.7321 

3.3798 

-2.2721 

-2.0447 

-2.0225 

1.7321 

69.3715 

-33.2071 

-16.1938 

-7.9688 

1.7321 
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/o-tl e °'' r 


anega = 
1.0e+03 * 


Columns 1 through 4 


2.1561 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Columns 5 through 8 


0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 


Columns 9 through 12 


0 
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0 

0 

0 
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0 

0 

0 

0 


0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0.2782 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 


0 

0 

0 

0 

0 
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0 
0 
0 
0 

0.1048 
0 
0 
0 

Column 13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 + O.OOOOi 
CXCLO-f- ^ 

X = 


/£-/? Co/j-f- 

0 

0 

0 

0.1805 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0.0500 

0 

0 

£Xaci Yf.f & 


0 

0 

0 

0 

0 

0 

0.0154 

0 

cf IS ' *1 


Columns 1 through 7 


10.1121 

0.0505 

12.3264 

0.1232 

19.9391 

0.2497 

36.2843 

0.4865 

68.5425 

0.9946 

132.3452 

5.9261 

358.7717 


-99.2702 
0.8057 
-60.8796 
0.9882 
24.5986 
0.4064 
91.0508 
-0.4901 
87.0711 
-1.0285 
15.2116 
-2.1565 
-103.5514 


-129.9937 

0.1529 

-123.0672 

0.2894 

-103.0258 

0.3951 

-72.0052 

0.4590 

-33.3062 

0.4871 

9.2843 

1.3595 

74.0764 


-75.3367 

1.2468 

-9.3374 

0.3091 

73.0222 

-1.1702 

27.4385 

-0.5987 

-66.2113 

1.0447 

-43.2898 

2.2411 

89.3228 


57.9828 

-1.2426 

-22.3957 

0.9599 

-40.6822 

0.5011 

53.8227 

-1.3465 

-0.8834 

0.5653 

-52.5323 

2.2287 

71.9851 


45.6155 

-0.8243 

-36.0458 

1.3027 

11.3519 

-1.2346 

18.1045 

0.6475 

-39.9837 

0.1777 

44.3715 

-2.1771 

-56.8220 


37.4875 
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-36.9790 

0.3350 
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39.1004 


Columns 8 through 13 
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-23.3039 

-19.2839 
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-1.1557. 

21.4342 

11.7057 

1.0915 

1.4029 

-16.1249 

5.0703 
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4ZD 


^ “ &er?ercL/ <f^ aai/on o£r>ioi?/On <o 

ISZ Ul 54 -/st o o"~] I 


u£ *£ L /*,( -st 1 o o j 

54- /}[ 212 0 S4 -ISt 

Z!L f L lt - _° 

0 0 54 ISt /56 -ZU 
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9, 


-.cm 

El 

->2 -It j M 0 -U 

j 


A 

~/.o 
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u u z \ o tt 1 -a zt L 


|^| 

0 


0 0 \~>2 ~ L l (iirjlQ-tt 


*3 

i 

-so 


o o 1 a *e 2 




■tessi 

L / 


W$t) 
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<n - 


*T\ 


Symmetric, er\ode (free-free V'i^rett/e/n) 

6 i = 0 e, = - Qj 

a re re ce$s a r y 


ue to 5 V rn me fry c »/- y Z~ & ^rd <tr 

/ ( / z 


~rom n be re ua t/or r>e hare 


/)*£ 

4ZO 


/£l 221 9t ' 

lit 4l l 13 £ 
£4 M 2IX 


C IT, ) 

e, 





o 

0 


Let *L = 2 
El 


ter soluc+io^ 4e>l(ouJS.' 


\*J 


|4U z z ‘M'C Sftr,y ><5 cxre •Aex iiJe and 4 lie. beam 

El / 

loeko.^9 albeit a. r^,J ^ jrUw/v, by 

4ke first m o Ae. 
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10-30 

cont- 


m = 

156 

22 

54 

22 

4 

13 

54 

13 

312 


k = 

14 6 -12 

6 4-6 

-12 -6 24 

»a=qln , *k*qiji 

a = 

0.0897 0.3402 -0.2032 

0.3402 2.8817 -1.3864 

-0.2032 -1.3864 0.7490 

»[v,d]=eig(a) 


0.8052 -0.1116 -0.5824 

0.1784 -0.8911 0.4173 

0.5656 0.4399 0.6976 

d = 

0.0224 0 0 

0 3.6087 0 

0 0 0.0893 


= iiv Ccho/ 0)1 J 
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z&jzij _ Ip of melton for B~x a/r)fJe. /O.S'—f /s 

4(2) £> <2 ) 


•— 


** 


3/2 0 

5* -f$t 

-X 

0 2£ X 

csl ~.lti V 



tSL -til 


-£.£“£ -nsi 

1 -t,l l l _ 


/'-T -^r *4 ~ 


%A- 0 —12 Si 

0 2i* ~2l >5t 

-12 Si *1 ~ 3 * 

2 l >£ l ~3t i. 



a/ here. / has been 
2 


a l^l \(JL ) = ? *?C 

A = W-fo /U ElJ <>724 El . r t 

W ' zuhst.tvted M . 

Sr nee i inf/e/e cf matrices f^ac/or oh./ of Me epuat/on 

Lei i = tO ">S/de of He matncet 


# 

Input ") to eo/upufer f>ro pra m 
Sc/uiion ep/veS 

&, - 3 • ^7 'f^£ a - 


(J u = 22 . 22. 
U> r 7S. / 6 

•J 


/- /. 


A-L, 


= 2Jg./ 
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/ 0 - 3 1 ^ orit 


m = 

312.0000 0 54.0000 -6. 000 

0 2.0000 6.5000 -0.7500 

54.0000 6.5000 156.0000 -11.0000 

-6.5000 -0.7500 -11.0000 1.0000 

k = 


24.0000 0 -12.0000 3.0000 

0 2.0000 -3.0000 0.5000 

-12.0000 -3.0000 12.0000 -3.0000 

3.0000 0.5000 -3.0000 1.0000 


a = 

0.0769 
0 

-0.0816 
0.3811 

»[v,d]=eig(a) 

v = 

-0.0542 0.1011 0.5904 -0.7989 

-0.2219 -0.8784 0.3807 0.1852 

0.2119 0.3722 0.7108 0.5580 

-0.9502 0.2823 0.0359 0.1267 


d = 

7.0810 000 

0 0.8406 0 0 

0 0 0.0018 0 

000 0.0735 


(jj ^ oil's'. Is 
7 - 3 ". / (c 
3. S~/ 7 



0 

1.0000 

-0.5996 

1.2864 


-0.0816 

-0.5996 

0.4580 

-1.3319 


0.3811 

1.2864 

-1.3319 

6.4619 
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rja o « -‘- s n 

. o It' -•»■** 

“'j " I 44 0 . 5 ? < 5 * " l| * 

[-t.« - nst -nt 1 

14 0 ““ “j 

4 - 0 U * "" I 

? 0 .$*£ ^ 


|kUU 

m is 


4 ,f»©*V 


<uu to L, CL. * * Q .; = * k " 4 1- 
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( 0-22 


>*n 


m = 

312.0000 0 54.0000 - 6.5000 

0 2.0000 6.5000 - 0.7500 

54.0000 6.5000 156.0000 - 11.0000 

- 6.5000 - 0.7500 - 11.0000 1.0000 

>>k 
k = 


25.8600 0.1310 - 11.3600 2.9220 

0.1310 2.0120 - 2.9230 0.4910 

- 11.3600 - 2.9230 13.8600 - 3.1310 

2.9220 0.4910 - 3.1310 1.0120 


a = 

0.0829 0.0052 - 0.0822 0.3817 

0.0052 1.0060 - 0.5986 1.2847 

- 0.0822 - 0.5986 0.4697 - 1.3314 

0.3817 1.2847 - 1.3314 6.4735 

»[v,d]=eig(a) 


v = 

- 0.0544 0.0949 

- 0.2215 - 0.8779 

0.2118 0.3746 

- 0.9503 0.2827 


d = 

7.0915 0 

0 0.8472 

0 0 

0 0 


0.6111 - 0.7839 
0.3740 0.2007 
0.6968 0.5739 
0.0331 0.1260 


2 A <4 

^ r iu ry\ H ^ 


0 

0 

0.0131 

0 


0 

0 

0 

0.0804 


do - Xl<€ 

is'. 7 r 
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0 ' 33 

>ra 


1.0000 

6.5000 

- 0.7500 

0 

0 


1.0000 

- 3.0000 

0.5000 

0 

0 


6.5500 

312.0000 

0 

54.0000 

- 6.5000 


- 3.0000 

25.8600 

0.1310 

- 11.3600 

2.9230 


- 0.7500 

0 

2.0000 

6.5000 

- 0.7500 


0.5000 

0.1310 

2.0120 

- 2.9230 

0.4910 


0 

54.0000 

6.5000 
156.0000 
- 11.0000 


0 

- 11.3600 

- 2.9230 

13.8600 

- 3.1310 


0 

- 6.5000 

- 0.7500 

- 11.0000 

1.0000 


0 

2.9230 

0.4910 

- 3.1310 

1.0120 


1.0000 

- 0.5822 

1.2269 

- 0.3532 

0.5713 

0.5822 

0.4016 

- 0.6448 

0.0946 

0.1294 

1.2269 

- 0.6448 

2.7725 

- 1.2088 

2.3238 

• 0.3532 

0.0946 

- 1.2088 

0.7127 

- 1.7975 

0.5713 

0.1294 

2.3238 

- 1.7975 

7.6275 


»[v,d]=eig(a) 


0.1311 

- 0.4785 

0.6513 

0.5105 

- 0.2625 

- 0.0268 

0.3343 

- 0.3957 

0.4519 

- 0.7258 

0.3888 

- 0.6818 

- 0.4418 

- 0.3241 

- 0.2893 

- 0.2458 

0.1319 

0.4262 

- 0.6481 

- 0.5661 

0.8778 

0.4206 

0.2057 

- 0.1004 

- 0.0133 


d = 


9.2416 

0 

0 

0 

0 


0 

2.7500 

0 

0 

0 


0 

0 

0.4707 

0 

0 


0 

0 

0 

0.0418 

0 


'y - rnJL^ 

07-30 Bj 

0 ^ 

0 / 25' (# 

0 sib. 5 j 

0 / (o / J- (o 
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Cjlobed . 
caorcLiwits 


Bj. 10.4-7 ~ K 
Etf. (64-£ — M 

B-le*nc.n{ ( 2 >~® ^ 

La Ad s jCL*** <ix 

*'"' W ^ * fafy/j*, 

* -22 ?••%/ 

i ./ < 

0 [ VC 

su = 2 **/ 1 y £'(* 

^' ~ d ~ Ii/[ *&'£*>■‘b 



o 'o 


-*-r f .:*-r ,w_^ 

Loa.4 S JTL 

VirtueU etrsf}. = ^ (X ) - £ U^t) 





Sujt ) 


r x I 

Su/= a-* 4 « jjtyc,; * -J p'fy/'U 1 

= 4 - [^*?;/£ ^ -y 2 ^ ■/*/£'**] 

-- ^‘4 - /*4 - 2 J 


See. Jzy, IO.Z-2. £or e.vcjucd/bn c-P / n feyr a. Js 


?<*) = 4 f (*> <r t i- 4 x (x)B ) + <p C x) a. r + ^ any *-lt.*><n't 

bed* */,t.) *nd x. 
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/o -M • 


>= m = [366 -22 0 0 ; -22 8 -13 
m = 

366 -22 0 0 

-22 8 -13 -3 

0 -13 156 22 

0 -3 22 4 

»k=[ 12 -6 0 0 ; -6 8 6 2 ; 

k = 

12-6 0 0 
-6 8 6 2 

0 6 12 6 

0 2 6 4 


=: ~fi~ “ o 

- 3 ; 0 -13 156 22 ; 0 -3 22 4 ] 

0 6 12 6 ; 0 2 6 4 ] 


0.0328 - 0.1068 - 0.0470 - 0.0711 
- 0.1068 1.0965 0.6856 1.3326 
- 0.0470 0.6856 0.4829 1.2102 
- 0.0711 1.3326 1.2102 5.4162 


»[v,d]=eig(a) 
v = 


0.5728 0.8135 0.0993 - 0.0176 

0.4220 - 0.1883 - 0.8417 0.2792 

- 0.7002 0.5464 - 0.3957 0.2340 

0.0602 - 0.0654 0.3537 0.9311 


0.0041 0 

0 0.0317 

0 0 

0 0 


Hq&I 

0 0 

0 0 

0.8714 0 

0 6.1213 


( jU 2 / ? 

3 . kVv" 
/?. /vj 
70 
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_f 0 ~ 3 ^ Conf 

Caeneral iz«4 -£©»*ee 


Element l"i 


— A * . A -£ 


T « S?7 -. 044ZSt 


~ 0 * 42 ?£ .otSSlO 
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e . 


Element 1-3 


Co / 


<msi.L 


I jV - * 1 [ # 

J • -. 0^*1 * 


must *■< mex \ y \ on 
vi ^ M - $\de 

force ^ moment 


c 

S u,(t>= -MtiL £ ^ j j <p.<^ i*. - -/Wlil/ 


To-feJ ^ti^Pnei^ 


I.X . l-£ — 1-2 • I i 

• 133 l* " **33* 
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it 
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Ob IX L 



(Vi ' 
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+ An£Lji 


*5 

r 

— 
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I 


t - 47 gr * -« 


•^SO -064if 
j-.&mtf . oi 3 ?j 




tr 

JL * 


~ ^ ** n * /'* 

L*vfu+ tnt® ct *® , j «-c . witU te "g-~ = I 


^ = 


Cj 


1 ml * 


AXD £l 
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/6- 3 V 


G o n ^ ■ 


r*\ 




w\ v _ f q 

bx. 




t/# d 


G as (2_ 


»k=[13.63 -6.164 0 0 
k = 

13.6300 -6.1640 

-6.1640 8.1570 

0 6.0000 

0 2.0000 

»a=qin' *k*qin 

a = 

0.0372 -0.1081 

-0.1081 1.1180 

-0.0476 0.6951 

-0.0720 1.3469 

»[v,d]=eig(a) 

v = 

-0.0178 0.0992 

0.2821 -0.8419 

0.2349 -0.3928 

0.9300 0.3565 


d = 

6.1429 0 

0 0.8846 

0 0 

0 0 


; -6.164 8.157 6 2 ; 


0 0 

6.0000 2.0000 

12.0000 6.0000 

6.0000 4.0000 


-0.0476 -0.0720 

0.6951 1.3469 

0.4871 1.2165 

1.2165 5.4257 


-0.5175 0.8497 

-0.4318 -0.1588 

0.7359 0.4989 

-0.0648 -0.0616 


0 0 

0 0 

0.0057 0 

0 0.0347 


0 6 12 6 ; 0 2 6 4] 


^ - Lti z rn 


cP ^ S~0- S' 

19. 27-(p 
3 .'S'if 3- 
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From E-q. (0.1-/ ArtA Sac. (0.4 "M*. *t<flfn«s£ 

for •fUe | ^ a\ A*m ent of Lar\ejt\y ^ 

• • 
at the pmv\£cl €.r\JL IS 

<Et “Mfo *(£l‘ 1M 

j£* “4(fc 12 ~4(1) OrA + K e, 

_*« o ‘ -*«) - ni)*J Ux . 


$a \(*tf> ' 2t '"‘If 6 ') 

' ~js -n 'z -*£ l 

I_.r^ ->< i x [a t j 

ro'to.'t tJL f 4-Uc «^er>er<a.|| ie et force <ft*er\ m £x.10.8-2 
£©r element «F i } requires ouAA\r\^ fkc Col \ 0 yu\r^ 

x -o*>/4(t)‘ - •*'+**&)-.MfJk+df (©/ ' 

** A 1 •♦*«• -°*u?(£) <"il 

(if Is , 
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• -\ _ 

A '/m>x3 x*7EI 

“H. ©fe>-t«.ir\ rciM.ce<i s+.ffne**, 566 ExiO.S" - * 

jS* rearrange rows w‘i+K a*?fkce"XSnb ,X -firs* Hires rows 
Ank. Angles .V> row S •4-, S", £ . 

Then rearrange ce>tu*r\ns k Conform t? neW column v£.c+*«* 
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r " aa O — 12. 6 00 

- 12. -4 **4- O _l2 - ^ 

_o_0_-IX ~~^> 11 ~ L 

~ o 0 -t X o O 

U X O * -(, Z 

o o i ^ -t ■4’ 


0 0 



ro r th<- twc 


reaurr 


*<L yy\CLiS m&rtnx 


?€«u-ran<j< Co l«.*rv nf 

*4.-11 0 | O 6 0 ^ > 

-i* x4* -a I -£ o £ *5 

o _— ( i_ ti I o —4 

0 -6 0 1 0 I ^ I 

4 0 — 6 ^ 1 £ ^ ■ 

0 6 -6l O 1 4 


Six 54 O ! o 

£4 -j i z £4- I $ 

_0 _ £4- j5"6i_0 

0 *3 0 1 g 

“‘3 o i 3 -$ 

0 -IS O 


-13 O K 

O -13 ^ 

_/3 JL2< ^ 

-3 0 

<0 -3 % 

—3 -4 


Kt V 


Assume. Mu — — Mi 2 — O 


M» _i_0 J V 

o TV "I* 


M„ V 


+- [JtLI* 

1 K21 ' Kzz 


\<^ ]/ = O — 3^3 4g q u , o y-i 


wWere «* = I<W ~ K j)L Ku Ki, = *+,£*«? 

Tsii 5"4 0 * 

Al'Hvow-^ M, v =^q£-| ^ J |4 p* wc . e*.n Choose discrete masses 

^ 0 ^4 I'S’t 45 £cf|ouJ$ 
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M. 


lfl-3% £gnt. 

ie. 055 U«e +kc «"*** of" CAck element At fheen^Li 


Ml 


~o r 


K 


K, t £ k „ 


133 ° 

© 

° ‘113 

o 

O & * 

l(,L>£ 

__ 

-1 

'*8 2 

X 

S 

it 

1 S 

1 

_0 3 

~ 0 6 0 1 

If* 

-(, O C, 

W-6 

n 6 — (o —lo 

JL 4 


.-f 

~xS 

-8 

4 ' 

= j- 

-8 

32 

-iU 

%0<g 

J j 

4 

-/4 

Cc> 


o 

2 

4 


-S 4 
32 44 
-IL <°0 


i r 


0 -L C 
L O -L 
° C, ~L 


_L 

Z.08 


-m xito 
~S’’j(. -17x8 



"s-.SJS ~I-‘^ AC ~*-‘ lL€ J 

JL = 

-S .308 


_ ii7 ^ -S.JeS fo. 3 ^ 



~%A 

-11 01 


'5*.^? -2.74f* 

•* 

-11 

X4- -111 

— 

-/.#*4 /3.844 - 8.? 08 


0 

“»*■ ix\ 


[_-1.7t1 10.3^ 



1 

fl«-46Z -10.4/S- 2.74-f 


-JO-M-T 10./54 


* 7ET 


*\ _ (eo* i*n£) _A_ 

A “ A7fcl 


o* inA” 4 
77 EX 


c0 


- 
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}b~ 3F + 

»m= [.333 0 0 ; 0 .333 0 ; 0 0 .1665] 
m = 

0.3330 0 0 

0 0.3330 0 

0 0 0.1665 


reduced 
(X> \ l-A. c\ j's 0 r-e,} q 


»k=[ 18.46 -10.61 2.769; -10.61 10.15 -3.692; 2.769 -3.692 1.615] 


k = 


18.4600 

-10.6100 

2.7690 


-10.6100 

10.1500 

-3.6920 


2.7690 

-3.6920 

1.6150 


4 




CD g> - 




% 




**4 


a = 

55.4354 -31.8619 11.7596 

-31.8619 30.4805 -15.6795 

11.7596 -15.6795 9.6997 

»[v,d]=eig(a) 


v = 


0.7873 -0.5901 0.1786 

-0.5634 -0.5710 0.5971 

0.2503 0.5708 0.7820 


oLT" G- L 


d = 


81.9738 0 0 

0 13.2281 0 

0 0 0.4137 


LP - 


V7-0 S' 

*.r? < ^ 
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JO'S# £0*1-. 


m = 


312 0 

0 8 

54 13 

-13 -3 

0 0 

0 0 


54 -13 

13 -3 

312 0 

0 8 

54 13 

-13 -3 


0 0 

0 0 

54 -13 

13 -3 

156 -22 

-22 4 


k = 


■u 

l£_ 


© 


C 2 > 


=J 

7 

/ 

X 


% 


% 





24 

0 

-12 

6 

0 

0 

0 

8 

-6 

2 

0 

0 

-12 

-6 

24 

0 

-12 

6 

6 

2 

0 

8 

-6 

2 

0 

0 

-12 

-6 

12 

-6 

0 

0 

6 

2 

-6 

4 


0.0769 

0 

-0.0545 

0.1691 

-0.0555 

0.0912 

0 

1.0000 

-0.4004 

0.7812 

-0.2136 

0.3415 

0.0545 

-0.4004 

0.2469 

-0.3327 

0.0045 

0.2650 

0.1691 

0.7812 

-0.3327 

1.9215 

-0.9207 

1.8330 

0.0555 

-0.2136 

0.0045 

-0.9207 

0.5897 

-1.5889 

0.0912 

0.3415 

0.2650 

1.8330 

-1.5889 

7.1247 


»[v,d]=eig(a) 


0.0183 

-0.0567 

-0.1449 

-0.6481 

0.7052 

-0.2409 

0.0837 

-0.5490 

0.7351 

0.2340 

0.2639 

-0.1636 

0.0130 

0.2971 

-0.2094 

0.5674 

0.2667 

-0.6889 

0.3126 

-0.6894 

-0.4922 

-0.0446 

-0.3072 

-0.2974 

-0.2321 

0.1647 

0.3399 

-0.4345 

-0.5128 

-0.5930 

0.9170 

0.3237 

0.1926 

2=- 

-0.1113 -0.0671 

co 2 -'* j2.^ 

-0.0192 

8.1884 

0 

0 

0 

0 

0 

0 

2.0602 

0 

0 

0 

0 

0 

0 

0.5817 

0 

0 

0 

0 

0 

0 

0.1147 

0 

0 

0 

0 

0 

0 

0.0144 

0 

0 

0 

0 

0 

0 

0.0004 


/V o. 

. iJ 

3 -S'/ cc 
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10 - hO 


!e 


*h 


JA 


4_ 

/ 

/ 

/ 




* 

3 


5 


/=/ 


Tht wass d/tJ slftf/res* mtrices \*r the elerrcsrh f# 

the fi'rsh ( krpr i> t/f'd^eh^) %tch* if ire. jfivefj cli 

j. f * n .j j£ r i 


2 

1 l 


Fo 


the 





■ L-l I. 

[Iifr stch’*», they ore y've/t fl-r 

qmJ . 8 . ft ~l 

L-i i. 


u,=o 

the 


7W. 

syste/t o.r-i 


art the 



1 0 o 0 a a o' 
$ 1 0 0 0 0 0 
x s 1 1 0 0 o 

o 2. * * 0 <0 0 

e o i hi o o 
o a o i ^ i o 

° 9 0 0 I k I 

0 o 0 0 0 I 2 


Mss 0/td S k''ff/tcss f»r~ 

JiVe/r as (Mfkck U{) 

r/fjO o o o o o 
-i If -2 0 0 o O 0 

Ft If ~XO 0 0 0 
0 0 -l 3 -I 0 0 0 

0 0 0 -I 1 -I 0 0 

0 o 0 o -» 2. -I 0 

0 0 0 0 o-t*~ 4 

0 0 0 0 o 0 -1 | J 


> 4 


Let’ U- [Hz fifj 11$ \lj Uf tig Ujl 7 . Therefor^^ 
the e^,UoJ:f'«r\ f t r the. knift uJtex! vtira-lhri /j 
Hii f ktl = 0 . 

Th t /Httur^l ftrcit/e/lctfr urt the square roots of tht 

yneroiiiieJf fr'Jtnyctlucs of ft a* J k J c-e.y k - =z 0 . 

S ol]//'/rj for' yj j Jtk 

uj — \tc, Ff. flSjIf ■ OlfOy IF [oiy 3. olfy/tafly f. 55ify, 3-130 

\i a noA~ 'Mve_ cursoo-t <-» &~d/y\ He 

/<**• Cua *«• , vclU?S 

proird^ F 0r ^Ke .-^ 00 reL ^.// 

Values 

P = ’C a / 2 / a. / 3^ £ O ^ "f^ 5 
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JO'*-/1 Co/vT oixvpor p;4~yea»i.y>) 

anega = 

54.5345 0 0 0 <2 “ e,l€mes>t~ 

0 18.7893 0 0 

0 0 0.8794 0 

000 5.5554 


anega = 


Columns 1 through 7 


538.6989 


0 


0 


0 


0 

0 

0 

0 

390 

.6837 


0 


0 


0 

0 

0 

0 


0 

219 

.6134 


0 


0 

0 

0 

0 


0 


0 

297 

.0578 


0 

0 

0 

0 


0 


0 


0 

160 

.7123 

0 

0 

0 


0 


0 


0 


0 

117.0056 

0 

0 


0 


0 


0 


0 

0 

75.8829 

0 


0 


0 


0 


0 

0 

0 

0 


0 


0 


0 


0 

0 

0 

0 


0 


0 


0 


0 

0 

0 

0 


0 


0 


0 


0 

0 

0 

0 


0 


0 


0 


0 

0 

0 

Columns 8 

through 12 









0 


0 


0 


0 


0 



0 

0 


0 

0 


0 

0 


0 

0 


0 

0 

G - 

g,le*ne' 

0 


0 


0 


0 


0 


]Q^.a.ry\ 

0 


0 


0 


0 


0 



0 


0 


0 


0 


0 



0 


0 


0 


0 


0 



50.7158 


0 


0 


0 


0 



0 

30. 

4202 


0 


0 


0 



0 


0 

15. 

4525 


0 


0 



0 


0 


0 

5. 

5100 


0 



0 


0 


0 


0 

0. 

8790 




<re.e /0 ~~ s,/ 
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10-12. 

... The mmi a/id sliffteu ^^{rTct J of 

Me ysh* art fll'vtft a $ (ij/hr* U, str/cl 

ft '- ?}u*l to Zero) let u= CU S J 




lAi. j- TlAz. 
A 

- -_L£l 

A 



2t.il/1 






El/KjA/ n i 


IT] 


U 4c 

J7 At 7 T 


T 


/Irina *t 

U , 


rr 


There fs also 
opposric t<> the 
111e t j Hahn>Y\ o f 


ci Jayi'/ij force at 3 e^Uaj f a ££/ 3 
*/'rcch*> of *ioho/t. There Are, 
motion fs y/rtihn 1 : 







For a st'jyle ']£= 1/ ~ 0 y u/Z Ioyc 

F,] = IlF u 


T i« *W/ %J ' 

ft =Q (/% i!fJ l ‘=3~> 9 l = ~F V > 

F, = JLl i 12 7 i f 4/<?,J = 3_£l jj (f CuJthy-r, 


J 


u/j 


Jflt HUM /j Pt>l/>1t] Ijf rtkll/jlv !/>j tfc P/Yjf r-JiAS 
of the »qj j W/-A i/ ^ *} 0, ^ ^ 7 r 

=2> Jdll l IS fiT/ t nJ §,] = M rtl nr 

HO L kO _ _ / 

Lf$ 1 Fi J' I ^l Y 

flunks g / fllolhY) 

1)41"/}- -mp -tyrfi) 

JtL pi-tf j-(F^4tj)y « fcy 

/>o /' i y 3 / 4 ' 



V 1 r * •' n i r r r r t< ~n~r-m TTT 


ItfO " ^ 7 / ' 

a) Jti i£z% 

; r w y 

=3 £ V£ +(k,t*c* W M 

3) a* J F + ( *>+f)?>= H 

I 

If y =y L = 7i kht forces fr 0 v I; ivl ^ ire?sr 0 . 

faree V/lr'oith>r fay Me khr?e l a i sm 3 
=^> The/ MJrt\ \/fir a.ie u/ fh tie :we fre^vp/icy UJ 


IM nf 
m 


C°yi I[ /or) , .! J/y■■’.»o / 7.C 7-Z_2i 

' t *y~ 


/y; 
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7 - / 


xk x .-y k x - K 

* ~ £ (I - c<rs cJt) , ^ = 


- - - 1? 
' ' ***** ^ £ 


f V ^ "/ f -- 5 

—£*■ cJ t - £ 

f* - 5 f' - -£=r -7 - *; j 


~tw*. 1L - 

/. y / - £* 


[/ - e 


/^/- /«*// C(«rr>f>,nj = Q 


-<*■ ~ , 


-r <j.^ 
<» 


C^5 oj 


t € j - f ^ 


. /w r' 


# ** 


f. ***«.•?.• *".■?.• - 

0 

where f(%t) - 3 f(t) Q.nd J> g - foto.1 force 

t. 

r t 

Molt pari/c/p a.t/on factor s -g ( <p. (.x) d%. 

O 


jLLz.. ,4 ^ ^ V ( / r ^ ) 

G*«*r*l.Z<J force = ^\j j roof to** } fet) 

(See Sec.//./ J 0 

H- * J J" -f>(x)^Xx) <*# , f>(x) - *■ S(x-e-) 

° r* 

•'• y l So-a.) 4>Xx)eLx - fj(a.) 

o 

?•■■ ¥t/ D -' rt; , -I ?,* - J iffafegfi, 


/*7 • 4 ; •* 

* 4 


* ct > 
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n-s 


K - l t~o iJ ]*<•>** 

O 

<f>. (u + e) - j _ £ d 4>" (X) | 


= Lm* t | 
£-rO L 


- p c • J~ d ^A* ) 
d? 


ch* 


X s CL 




= fci 4 >.(cl) 

Ste A r r*nd. D f° r <P 


LLzJb - 


*i = if t - taj 


<f> (v) - fT S~* «.ir £ *'>*'> normalized to total ">** M 

<j> (l) - /a" S*~ **• - / y S f • * - 

u)‘ = />„*)* ££. /t = *-? ^ " 

~ {P ~ ' nt* l 

<f,(i* t) - y ( 

7 ' & £ 1 


‘-'"•f-s—Zl DCV 


m.z / S i> 
• > / 


J/..--L 


Fr cm Kt. n-s x m ‘fe£'<v /.= f 


4>(X) - yr 5 Kir *■ - f£ if ^ *.7T£ 

~ dsK ' l d 


mir 


s % (o 


s'/ 


<d)‘t t Jij * r - . 

+“-f~ = f/»-< iVK*; 

r- * 3k? ^ *'0^ 

* ( ,, t ; . J ^ ■ §/l 7^ ^ 


M 

x _ c< 5 JS2* W 

-Z —^ 

m.r a • ' • 
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//- ? 


C)n/y r modes #1—2 £ - wtf/ atre. 

/ / ' * J 


k £ o 

m, ' 


ji =)Ts***^ 

-<4 


* =f 


0 ^ 

P * « J"cJ, (t-f)ctj = /- 

o } 

?(',». z *4^ ^ -L^ 

m* 2 C ’ - - 

/ * 

- LhL y Z (/- c-3 ST) W "^2 £ A - OJW.t ) 

-7T/V7 Z_ a- / ^ 1 ~ ' 


*n = 1 c 


, 5Z 1 


/ »r ode s i 


« (* t) = f s<^ *£* rI- cjo^-t) 

*' ' -TTMCJ*- l ' ' 

A 


mode, *>*(’ 


■ ££3. C'-^O 


//- i 


fi/ormul h*ode U (X) = G. ^ +l3 <*JT) *>* 

If -g{ <2e -* -*»*-*? ^ 


Fof") 


4 


~3our>d*ry C»n d. 

x = 0 da. ~ o 

' 4* 


a = o 




x - l , u. - q /, t3 - o 

u)t m ir jv s^T - - - 

C X , a. , a 


/j s ^ * £ fdO + c**^)** = ■£ 

o * 


• '. p - fT ^ «f| 


■n. = / S .T 

' > y 
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//- 4 : l 

- h ^ <re « rf * = ^ ' 


f- 


__ j, = /r^. b w = >£ rHt'-ZLV 

h C^M * ‘ (*££)*" 

u ( *t) s T& l *>t~w * —i c l5±J > ' < * ) + tli£ ^ (t) . 

' «■ 4 ^ l r?)*' dr*J x 


//- /o 


K 


•~ '- i j‘ s(*-f)rT'”¥ : z •** */*«» 


2l2T* 

£> 


~ / S &* 


ov 


/ 

* 

C<o JO * 




J 

Co ^ = 

C*> ?0* 

= O 

mode OL^Sent 


5 

= 

6 

Co /'^T0° 

- -.?<c 



7 

eo 7 — - 

W 6. 

d*D X.)0° 

- -.St4 



9 

IT « 

<# 

C*. 5 Z7D* 

- 0 

mode ^Stnt 


// 

^ «f - 

CC S3 o' 

• = .stt 



'/ O 

des present 

are. / 

y 

iT 7 // 

/ ' ' 

/3 - - - 

/ 


// 

- /! 

= 

|T Co **£ 

Vi 







f- 

- Fo 

U> m l M 


X 

^ =■< 

f *k7T C \* _ / 1 
i X. Z/ " V X > 

' M t 

m. C*, 

,t) - 

/4f 

1 (f/ 

Z£ 

_ Co €212 

(fO" “ 

■Pr «; 

- 

g-Cii iv 7 Jfr.2>(t) 
(TJTj* 7 

*• £Jii£ c«5 

III-TT ) x 

1 t ' 

• ?,«-> - ■ 

-j 


346 




Assume. = (£ m 


£ u 

U*[iii x +icrfifi]t T = u f'“ . 


.1 ,.4 

sr x. 


<j*, ±. b L -+ i *V 
Moment*-* o S>'*S ~ x ~ ^ j-j 

^-«^^t4j* 80 .*. U -m<S-(-zL mbi) = 


* *£ - *'/*, 


«9/- tO ' X 


r # " 7 / 

<*U*S*»c ,*•/ ^ , *■ /A . a. 

-#4- = ■:V>y- /. i\?o-x)j>+ttx =o 

iS -t+b +4 ^ 

7 r *?- 


^ = £[(*-0 ± f 7 *~ 0 '+ & r ] 
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- /J 


r2t?. 


Cfc 


£cntrc^//3.ecl £orze. 

,-e 

= %f(t) / -£)<** (i)fto 


id 9-2. Hl)U‘> u " e ?•*" ° r 

* ' / H ° l T* , M S-> 


Sol ution 


^/Vv 

(“rom Ar>pCrv^ix"D ° 


4>(£) r /,5“SS ^ -4.73 4>,V> - X.O 

<£({) c O fct* 7 -« <*■>> = J ° 

i)>(i ) = -|. 371 fyt * I0.*n ^'(OJ =1.0 

J 

'I. = = Hi 


11- !_± 



K- = 2 f ** 


The CthoS*. Jnie^ru.1 Tor f<£^ Can Le e xa ,/<* * bed Tro*r> “ta.bJes ' *> 

App. T). Add co f u. m n of ) untf rn u.l{„ jby /n?er*rod d£=.0.CP4 

/-mode Co/umn 5 <-"» /* ZD. £7 /, O-^ - 0, <?,Z;2_ £x+.ct 

vcJoLt from tnfcjrci.tinj £cf t c f^^I jL w//h yroyer boundary Candy t/onS 
« o.m ■ <*■** *(*) zCctiJi!,,.-<>*>**)-{ 


vcJfxt from / nfeyrc-tt nj £</,*? ^-/Z. w/Vh yrayer boundary Candy t/onS 

,s o in ■ '■'■“* * <*> * c**b - %-**-**) 

Mptt: The number ZD. £7 ts tv tar^e $/hce it u/clS corny u/e U With the 
tarqttf def/ectron of 2 . 0 x . o* + /.&?x • 04 •+ - * • etc. If the fur* >s 
Started *Jith /.^x^O^ + ••-etc ^ then a f/yuro of /&. £J fS 
obte^td. A^trayncj £ (/$ .£7 1 - *. 0 .£ 7 ) x * 04- u/e obtain the Cx*-*t 
va/ue 0.125. • , _ _ 


'betfiny o(^s — 


\l-S»Zi 


CozJLfl* c<Ap£ , 


the inte^rc^l Can be shou^n bo 


k*) & = J4" 

K 7 fa? 
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//**/*» */• = -*t ( 4 ) 

«(*,*) *«($, f; 1 = m h 

Fro** Etj./l'^ — ^ 

t(f) * -wt* rr +(eSi<0> i- ftyy*;* •• ‘ 

U$*nj J^mocfe o*/y 

«i)i> g] - m'%, *,d) 

pfj) - id )■*■ A, V^ r rAod*S <^A/c/ <rfi hej/fcictf. 

> C*~*r« rt *°* / -> 

t /co ^ r 


(& \ _ , + -^- .££* 


IS. = 0.0 /£~+ 

17“* 

±. =0.0/ 
x+3 
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//- /g 


Fr 0 ~t E<fj/.2-4- 


4 , ro = *(*, *-) ~~ (cJ/) f'^ f/ 0 -) 

*,<*>%'“’ = '^ w ?' w -( ! 8 l )V lt> *> < * > 

■■ (sS = ' *- A. ^ 


LL - a- 


Firey t Ctj, ieco 


f,r^-*•>*;= - £ fa) {%.?, '(*> + ^ <p>;J 

ojuj:^) = -* ^(°)[%, fa) + f *. ?j°) j 


% fa> ] 


M fa) fa) 


% fa>% '( 0 ) [ c fa fa + M fa)} 


- o 


\ - /c. oj* F< 


i t 

5/~f/>/^, Su.pporJzU beam t /.et \-(%J f ) 

[u -»+' w *sk>J 1 s ° 

j 1 

A a '-{ n + W * *?V>J J A i- {it *'4 *pM j “ ° 

f^») = fZ^ 7 ^ fa) = {l £ ?/%)**(& 

%(*) ~ f* ^ z j = a/T fa(<=>)~ ?(l) 


If X = © 


^ - / 7 A +■ IL 


ff X *0 L*i Oi = ^ • 


-^-/XA-/i) - ° 

A/ = ' 
Ax * tC 


f 


^ — [/7i-/o(-f : Jo<;J> l + [/ir 4o (^) ^ ] ” O 
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LL-Zi — C.oni: 


M. ‘ t 
0 


t - + ^ + [**?%!’** +i *o % jc 

0 * o 4 

= * f. f«,r ■*■ ,? * 3 } * **’* ?/* 

U - if. lei y'cxt)d>c - i f z h ^ * f/ 2 '| 5*« //*/- 7 

•© 

LoijrQ n^CS £< f*- 


Fre -v £c t 


%o(‘ Wo+i’i^J + x.*.12S*«l %/ — ° 

p lm e t- - 0 


— ^A 1 0 ^ 

-X> «<«' 


-x*.i r gj 

Z (uf-i-f)**!- 


- o 


cJ _ / ^c,/W r *7 t_ K/A<re U>, - / * nat. fr«y. of Canf/e. 

ZJ, ~ Jm /W-f .3?7 M a of l*nj i and *r)A& £ *7 


If to.-* ° , /Aen 
5 ,„c€ u>, = J.sz~f£f 


U> - CJ. 


- l.Lj uz. 


' fsti “ "•' 

<*„</ I.r free- free £&*"> of Le*jM\ -Z-£ 

A*s W. A# a; - ZXt+f^£— - r rr, fW 
J 'ff i 


(zef * s ' S 7 /„s~ 


£ltl s. J• EY ' Sioce Me ct irpJan£ bee.mts °- f r e*~free be*'* 
3 " S ~ X ,f? CenfM 3Lt *f to).-* 0 , af^ve resuJf ,s 
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LLzJA lonif. 


U » if 5 " £ r // ’ / ‘ 7 


L >4 j rc 4 * ^ es JEq ^ 

i 


C -4 a rcLc/e rt$ tt c 

k ** 






'x.ffdv 


'C/ c 

X'm(lcS‘ - cj-j y?^ J 


- 0 


/"o'* ca* 




- ^ 


z - z4 


*% 4 x- - ~ J iS - /C /. 7 S-? 


cf r 


\„e(L _ 4=^>‘ j = < jV W (ft */”'J 


cJ - aJ 


7 lfc t „t>)-£££L 


'*c 4 


me 
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0.2800 

- 0.0008 



0.3000 

- 0.0009 



0.3200 

- 0.0010 



0.3400 

- 0.0010 



0.3600 

- 0.0010 



0.3800 

- 0.0010 



0.4000 

- 0.0010 



0.4200 

- 0.0009 



0.4400 

- 0.0007 



0.4600 

- 0.0005 



0.4800 

- 0.0003 



0.5000 

- 0.0001 , 



0.5200 

O . OOOl ^ r - 



0.5400 

0.0004 



0.5600 

0.0007 



0.5800 

0.0009 



0.6000 

0.0011 



0.6200 

0.0012 



0.6400 

0.0012 



0.6600 

0.0011 



0.6800 

0.0008 



0.7000 

0.0004 rrX 1 



0.7200 

- 0.0004 



0.7400 

- 0.0014 



0.7600 

- 0.0028 



0.7800 

- 0.0046 



0.8000 

- 0.0070 
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^This is the Holzer method for the sping itbss sysyton 
clear; clg; 


m=[l 1 1]; 
k=[l 11]; 
for 1=1:13 
x(l)=l; 

w(I)=(1-1)*.2; 
larrteda(I) = (w(I) ) /s 2; 
x(2)=l- larribda(I)*m(l)/k(l); 
f (l)=larrbda(I)*m(l) *x(l) ; 
ff=f(1); 
term=l; 
for ii=2:3 
tem=tenrvfx (ii); 
f (ii)=larribda(I) *term; 
ff=ff+f(ii); 
x(ii+1)=x(ii)-ff/k(ii); 
end 

x4(I)=x(4) ; 
and 

plot (w,x4) 


[w 1 ,x4'] 
ans = 

0 1.0000 

0.2000 0.6111, 

0.4000 -0.4249^' 

0.6000 -1.7395 

0.8000 -2.7949 

1.0000 -3.0000 

1.2000 -1.8708 

1.4000 0.7617^^*- 

1.6000 4.4980 

1.8000 8.0710 

2.0000 9.0000 

2.2000 3.1993, . 

2.4000 -15.4598 3 


prob. 12-24 
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/2 'JS 


/R-a.7 

(j^e 'K. (f Sa.'KC. CX.S 

€LYhe.f\d ci i ndtx A) 


prob 2-27 



omega 


[w‘ f4 • ] 
ans = 

0 0 
0.2000 0.1365 

0.4000 0.2881 yj 

0.6000 -0.1475 1 

0.8000 -1.6546 
1.0000 -4.0000 

1.2000 -5.7559 

1.4000 -4.3056 

1.6000 3.2041^ t ^ 2 - 

1.8000 17.3942 

2.0000 32.0000 

2.2000 26.2565^ 

2.4000 -46.553:T ^ 




'J-t/. nn ' 


i3l- 3l<? 


Use 'M.c- pryr** 

5“ - Crc # /J*, '*r, /«*/ 'A{/yooo<$j 
/< - >0**C/S~ •?©, A 2^ •?<*, /io, 6 7 

TXje s^ste^i w c^s 0 ui ai 0 enitr /.',?€. u.s/^^. */ x 0 -£ 

'f'A.g, f-ixse ^ a.^ £- • ^/W-r-sy./7i /»? e./r'T c. ^ .s &/*■€. 

fi.)l/e>\ Po' ^(p = 'Q’<SLnd 5^ m r7ltfa~> <L, /7l ocA f\j p0/~ 


T 


0 . A'd^-e: a^e dlos*. P-tr 


C. &.cU O^er ^eC-G.us'C PA_g. i' < C '•J ^ 


»[w‘ theta6’ t6‘] 
ans = 

1.0e+10 * 


0 

0.00000000100000 

0.00000000200000 

0.00000000300000 

0.00000000400000 

0.00000000500000 

0.00000000600000 

0.00000000700000 

0.00000000800000 

0.00000000900000 

0.00000001000000 

0.00000001100000 

0.00000001200000 

0.00000001300000 

0.00000001400000 

0.00000001500000 

0.00000001600000 

0.00000001700000 

0.00000001800000 

0.00000001900000 

0.00000002000000 

0.00000002100000 

0.00000002200000 

0.00000002300000 

0.00000002400000 


0.00000000010000 

0.00000000009955 

0.00000000009819 

0.00000000009594 

0.00000000009283 

0.00000000008889 

0.00000000008416 

0.00000000007870 

0.00000000007256 

0.00000000006581 

0.00000000005852 

0.00000000005078 

0.00000000004267 

0.00000000003428 

0.00000000002570 

0.00000000001702 

0.00000000000834 

-0.00000000000023 

-0.00000000000862 

-0.00000000001672 

-0.00000000002446 

-0.00000000003173 

-0.00000000003846 

-0.00000000004459 

-0.00000000005004 


0 

0.00070454393355 

0.00277996177187 

0.00611264965261 

0.01051670778083 

0.01573903028613 

0.02146629176526 

0.02733369347504 

0.03293529614420 

0.03783573271471 

0.04158306348715 

0.04372250859723 

0.04381076893278 

0.04143062692522 

0.03620550350001 

0.02781363719313 

0.01600154634395 

0.00059643562529 

-0.01848278581237 

-0.04121619545193 

-0.06747556083990 

-0.09701933458121 

-0.12949019031575 

-0.16441531332092 

-0.20120961780241 
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LO 

0.00000002500000 

0.00000002600000 

0.00000002700000 

0.00000002800000 

0.00000002900000 

0.00000003000000 

0.00000003100000 

0.00000003200000 

0.00000003300000 

0.00000003400000 

0.00000003500000 

0.00000003600000 

0.00000003700000 

0.00000003800000 

0.00000003900000 

0.00000004000000 

0.00000004100000 

0.00000004200000 

0.00000004300000 

0.00000004400000 

0.00000004500000 

0.00000004600000 

0.00000004700000 

0.00000004800000 

0.00000004900000 

0.00000005000000 

0.00000005100000 

0.00000005200000 

0.00000005300000 

0.00000005400000 

0.00000005500000 

0.00000005600000 

0.00000005700000 

0.00000005800000 

0.00000005900000 

0.00000006000000 

0.00000006100000 

0.00000006200000 

0.00000006300000 

0.00000006400000 

0.00000006500000 

0.00000006600000 

0.00000006700000 

0.00000006800000 

0.00000006900000 


-0.00000000005475 
. -0.00000000005869 
-0.00000000006180 
-0.00000000006406 
-0.00000000006546 
-0.00000000006598 
-0.00000000006564 
-0.00000000006445 
-0.00000000006244 
-0.00000000005966 
-0.00000000005615 
-0.00000000005197 
-0.00000000004721 
-0.00000000004195 

-0.00000000003628 

-0.00000000003030 

-0.00000000002411 

-0.00000000001783 

-0.00000000001157 

-0.00000000000545 

0.00000000000042 

0.00000000000592 

0.00000000001094 

0.00000000001539 

0.00000000001916 

0.00000000002218 

0.00000000002436 

0.00000000002566 

0.00000000002603 

0.00000000002546 

0.00000000002394 

0.00000000002149 

0.00000000001817 

0.00000000001403 

0.00000000000918 

0.00000000000372 

-0.00000000000219 

-0.00000000000840 

-0.00000000001472 

-0.00000000002095 

-0.00000000002688 

-0.00000000003227 

-0.00000000003690 

-0.00000000004056 

-0.00000000004302 


-0.23918201574982 
-0.27754480972827 
-0.31542622451125 
-0.35188603048073 
-0.38593414575746 
-0.41655203470013 
-0.44271664844260 
-0.46342657933742 
-0.47773002645629 
-0.48475409468513 
-0.48373487657210 
-0.47404769518270 
-0.45523681914657 
-0.427043 « QCn ^ 
-0.38943432166274 
-0.34262062466385 
-0.28708209408002 
-0.22357962439411 
-0.15316492796073 
-0.07718318108356 
0.00273177604275 
0.08467051902641 
0.16647015791455 
0.24574615238018 
0.31993434882400 
0.38634359875886 
0.44221911472193 
0.48481647616685 
0.51148590960207 
0.51976613076509 
0.50748664777117 
0.47287697868056 
0.41468073033284 
0.33227191294532 
0.22577022201726 
0.09615130046836 
-0.05465280558446 
-0.22366866991186 
-0.40681863857555 
-0.59888783801372 
-0.79352994513813 
-0.98332191613568 
-1.15987927799034 
-1.31404511781761 
-1.43616756696207 
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Jol - oLO • 


UJ 

0.00000007000000 

0.00000007100000 

0.00000007200000 

0.00000007300000 

0.00000007400000 

0.00000007500000 

0.00000007600000 

0.00000007700000 

0.00000007800000 

0.00000007900000 

0.00000008000000 

0.00000008100000 

0.00000008200000 

0.00000008300000 

0.00000008400000 

0.00000008500000 


-0.00000000004411 

-0.00000000004366 

-0.00000000004159 

-0.00000000003784 

-0.00000000003247 

-0.00000000002563 

-0.00000000001759 

-0.00000000000878 

0.00000000000019 

0.00000000000849 

0.00000000001506 

0.00000000001853 

0.00000000001722 

0.00000000000905 

-0.00000000000848 

-0.00000000003835 


-1.51648237671449 

-1.54561912796732 

-1.51525171589879 

-1.41891601085994 

-1.25302002629042 

-1.01807453222036 

-0.72017484749937 

-0.37276753436557 

0.00126108536749 

0.36713426279700 

0.67532803833047 

0.85826050744362 

0.82647518628648 

0.46426215681840 

-0.37534396649941 

-1.87625512220048 


-1Z.-3Q 

l~ro~t con SCrvoiton 
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1 
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I 


%This pogram uses M/klestad's method on a three mass beam. It is 
% set up for evenly spaced lurrped masses, 
clear; clg; 
n=3; 

mass=[l 1 1] ; 

1=1; 

ei=.le7; 
coef=l/ei; 
coef2=(l~2)/(2*ei); 
coef3=(l /s 3) / (3*ei) ; 
for j=l:200 
anega(j)=150 +10*j; 
va(l)=0; 
ma(l)=0; 
thetaa(1)=0; 
ya(l)=l. ; 
for i=l:n 

va(i+l)=va(i)-mass(i) *amega(j) /s 2*ya(i) ; 
ma(i+l)=ma(i)-va(i+l)*1; 

thetaa (i+1) =thetaa (i) +ira (i+1) *coef+va (i+1) *coef 2 ; 
ya (i+1) =ya (i) +l*thetaa (i) +rra (i+1) *coef 2+va (i+1) *coef 3 ; 
and 

vb(l)=0; 
rrto(l)=0; 
thetab(l)=l; 
yb(l)=0; 
for i=l:n 

vb(i+l)='do(i)-rrass(i)*cm^a(j) /v 2*yb(i) ; 
rtb(i+l)=mb(i)-'do(i+l) *1; 

thetab (i+1) =thetab (i) +mb (i+1) *coef+vb (i+1) *coef 2 ; 
yb (i+1) =yb (i) +l*thetab (i) +rrib (i+1) *coef 2+vb (i+1) *coef 3 ; 
and 

thetal=-thetaa(4) / thetab(4) ; 
y4 (j) =ya(4) +yb(4)*thetal; 
and 

[anega 1 y4‘] 
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I3l~3) C«^h 


Q/hea-o- 
l.Oe+ro * 

P 



0.1600 

0.0006 

0KC<f4 

1.7200 

-$.$005 

0.1700 

0.0006 

1.7300 

-0.0005 

0.1800 

0.0006 

1.7400 

-0.0005 

0.1900 

0.0005 

1.7500 

-0.0005 

0.2000 

0.0005 

1.7600 

-0.0004 

0.2100 

0.0004 

1.7700 

-0.0004 

0.2200 

0.0004 

1.7800 

-0.0004 

0.2300 

0.0003 

1.7900 

-0.0004 

0.2400 

0.0003 

1.8000 

-0.0003 

0.2500 

0.0002 

1.8100 

-0.0003 

0.2600 

0.0002 

1.8200 

-0.0003 

0.2700 

0.0001 

1.8300 

-0.0002 

0.2800 

0.0001 

1.8400 

-0.0002 

0.2900 

0.000<£_, 4 JO 

1.8500 

-0.0002 

0.3000 

-0.0000 

1.8600 

-0.0002 

0.3100 

-0.0001 

1.8700 

-0.0001 

0.3200 

-0.0001 

1.8800 

-0.0001 

0.3300 

-0.0002 

1.8900 

-0.0001 

0.3400 

-0.0003 

1.9000 

-0.0000 

0.3500 

-0.0003 

1.9100 

-0.0000 

0.3600 

-0.0004 

1.9200 

0.0000 

0.3700 

-0.0004 

1.9300 

0.0000 

0.3800 

-0.0005 

1.9400 

0.0001 

0.3900 

-0.0005 

1.9500 

0.0001 

0.4000 

-0.0006 

1.9600 

0.0001 

0.4100 

-0.0006 

1.9700 

0.0002 

Vx£> \ ^' 4200 

-0.0007 

1.9800 

0.0002 

^ VeV ^ 0 - 4300 

-0.0007 



t 0.4400 

-0.0007 
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12-32 
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12 -il 



E^S. 12-4 ~ / % (o Can be drranjeet in 
■fol/ov'/nc yyi 6 .tr/jc /a c~m 
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J 0 
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/ 
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0 

V 0 ^ 


%) 



deiermt w*nfc 
to Z<.ro 
nat. £ . 
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coa^- 


%This program is for problem 12-36 and 12-35 
clear, clg; 
i=sqrt(-1); 
m=4; 

L=60; 

j=[20 100 10 30 ]; 
k=le6*[0 0 5 10]; 
c=[0 0 0 0 ]; 
g=[0 1 0 0 ]*le4; 
for ii=l:L 
w(ii)=ii*10; 
for j j =1:4 
zr=-w(iir2*j (jj); 
zi^w(ii)*c(jj); 
crj(jj)=zr+i*zi; 
zr=k(j j); 
zi^w(ii)*g(jj); 
crk(jj)=zr+i*zi; 
oil 

de(ii,l)=l; 

t(ii,l)=crj(l)*de(ii,l); 
for jj=2:m 

de(ii, jj)=de(ii,jj-l)+t(ii,jj-1)/crk(jj) ; 

t (ii, jj)=crj (jj) *de(ii, jj-l) + (l+crj (jj)/crk(jj) )*t(ii, jj-1) ; 

ad 

de4(ii)=de(ii,4); 
t4(ii)=t(ii,4); 
ail 
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/2.-V' C^- 

%This code is for problem 12-41 

clear,clg; 

m=2000; 

cm=50; 

st^> =20; 

j=t.2 .2 .2 .2 12.3 .168]; 
k=le6*[6.1 6.1 6.1 4.5 .0017]; 
for ii=l:m 
tl(ii,l)=0; 
theta(ii,l)=l; 
onsq=arT2; 

tr(ii,l)=-cmsg*j(1); 
for kk=l:5 
L=kk+1 

theta (ii, L) =theta (ii, kk) +tr (ii, kk) /k (kk) ; 
tl (ii,L)=tr (ii,kk) ; 

tr (ii,L)=tl (ii,L) -cmsq*j (L) *theta(ii,L) ; 
end 

cm=crrH-st^); 

and 
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We n«e.<X only 4o ctJc-.le.-te 4 he lo-st- two col^-nns 
O f' 4 |a€. (.&.$■£ tvJo v~owwfJ y w»4>cli CL re ,* 
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13 -/ _ 

(/j noise /« radio rece.pi.ion — •*//</<* frey. ranje ** />r*Aa bJy 

S tat/ on ary . 

(Z) ground rnof/on of ea r/AjuaAe — nan s/a/yen ary 

(3) ocean K/are be/fh/f — nor StaJ/onary and deyCnate y/t 

on m// n O'* S€o. 

/3 ~ 2 See text section 13.2. 


/3 3 See text sections 13.2-13.3 for ctefcn/Y'fin 

£(&) f°** $00 Co/ns ■/Arvu/n /cro d"»es 
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0 f zero or tnf/Ji'fe. f*.o.Xs *■* ■ z - e - ro 


/3-g 


4 r^orrya.1 pro bo.* 


,/.Sy 


e.ur t/e. 


= /.<? 


* 

i<r* 


<r fvir 


dpc = 


cr~ /g - 
(F fifr 


- /. o 


Ca.J 


405 



J3-9 CW - 
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7 * <r* 

•° 1 

f - <Z oipc ------ (bj 
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by /ntpvc/f* h. (& ) ~ O /or ^7 oe/ct. 
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*(/ + '/)■=■ ACat(t*t) 

c £ co>/ cc 7" - •£**» ^ ^J 


%(t) X(i*T) - A {ce^t Co 7"- Coy t SA./ J 

fy f Colt) - SA«y • S»A /to/J dt 


- r /z 


= /~ £[<*>• nu-*“•'>■■ wJ 


£0; 


T-+*o 
. i 


= A c„r 

X 



407 



J- 3 ~ Th/S proiJe~< /s Sm/Jor to Trol. /3 - # 7»/*jra te tAe area. 

under pet*?) curve , Start c.urve U/J-/A &CT) =. S at 

0-= O ctnd Lffear/y cte.cre.ase to &Cr) a f.o at 'T •= / etc.. 

Shift traced Curve as Sujjes/eU. 
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^ D J'O >2.5* 
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p- ='fet«4 Area = 1424 RMS = 37.9 */£*■ 
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/ 3 ~ ^ - ~p 
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/3-Zl CW, 
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In/tjr^-t'f Ar»<A.n\J conlour 


= XTTi V Rescues 
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-L^z32. Cj^ojL 

Ey c / f c^i / c n r m *n 3 (uj) 


sr«j)* £• + J = zl + W*£ \* s^f 

^ ^ Jr, ^ * Z/JJ 

&'-■)] 


= 10 


_/ 

L * 


jL + J- 

'rr'cj 3 - )\ *- •*- 


~ + O 1- _ 
Xo 


jportse 3f>ec i-r u m 
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zJ d£ 
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&=.to*(zrr) =37;i?J. ^ 
-Ht'-)43L. x /0 6 a/^,*• 


2 
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»tf=10 
tf = 


function xdot=probl441(t, x) 
xdot=[sin(x(2)) ;x(l) ] ; 


10 

»x0=[0, pi/3] 
xO = 


0 1.0472 

»[t,x]=ode23 ('probl441' ,tO,tf,xO) ; 
>5plot(t,x) 


prob. 14-41 




7 "V2- 


inHaJ CA 


as /*/i fj/'&ble*n 


function xdot=probl442 (t, x) 
xdot=[sin(x(2))-.3*x(l);x(l)]; 


prob 14-42 







clear 

global m g I FO w 

'Problem 14-44 Solution' 

m=input('Enter the value of the mass (kg) : ') 

l=input('Enter the value of the pendulum length (m) : ') 

FO=input('Enter the value of the forcing amplitude (N) : ') 

w=input('Enter the value of the forcing frequency (rad/s) : ') 

xO=input('Enter the two dimensional initial condition [x xd] : ') 

tO=input('Enter the value of the initial time (s) : ') 

Tf=input('Enter the value of the final time (s) : ') 

g=9.81; 

xpO=[xO w*tO]; 

T0=0; 

[T,xp]=ode45('p1444d',T0,Tf>xp0); 
m=size(xp);n=m(1,1); 

j=i; 

yi (i )=xp(i ,i); 
y2(1 )=xp(1,2); 
plot(y1(1),y2(1),'.') 
hold on 
for i=2:n 

xpo(3)=2*pi*j+xp0(3); 
if xp(i,3)-xpo(3)>0 & xp(i-1,3)-xpo(3)<0 

j=j +1 ! 

s1=(xp(i,1)-xp(i-1,1))/(xp(i,3)-xp(i-1,3)); 
yi (j)=s1 *(xpo(3)-xp(i-1,3))+xp(i-1,1); 
s2=(xp(i,2)-xp(i-1,2))/(xp(i,3)-xp(i-1,3)); 
y2(j)=s2*(xpo(3)-xp(i-1,3))+xp(i-1,2); 
plot(y1(j),y2(j),'.') 
hold on 
end 

if xp(i,3)==xpo(3) 

j=j+1; 

y1(j)=xp(i,1); 

y2(j)=xp(i,2); 

P lot(y1(j),y2(j),'.') 
hold on 
end 
end 

xlabel('theta');ylabel('thetadot');title('Problem 14-44:Poincare 


Section') 
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clear 

global m c k mu F w 
'Problem 14-45 Solution' 
m=input('Enter the value of the mass (kg) : ') 
c=input('Enter the value of the damping coefficient (N.m/s) : ') 
k=input('Enter the value of the spring stiffness (N/m): ') 
mu=input('Enter the value of mu (N/m A 3) with its sign : ') 
F=input('Enter the value of the forcing amplitude (N) : ') 
w=input('Enter the value of the forcing frequency (rad/s) : ') 
xO=input('Enter the two dimensional initial condition [x xd] : ') 
tO=input('Enter the value of the initial time (s) : ') 
Tf=input('Enter the value of the final time (s) : ') 
xpO=[xO w*tO]; 

T0=0; 

[T,xp]=ode45('p1445d',T0,Tf,xp0); 

m=size(xp);n=m(1,1); 

j=i; 

yi (1 )=xp(i ,1); 
y2(1 )=xp(1,2); 
plot(y1 (1 ),y2(1 
hold on 
for i=2:n 

xpo(3)=2*pi*j+xp0(3); 
if xp(i,3)-xpo(3)>0 & xp(i-1,3)-xpo(3)<0 
j-j+1 i 

s1=(xp(i,1 )-xp(i-1,1))/(xp(i,3)-xp(i-1,3)); 
yi (j)=s1 *(xpo(3)-xp(i-1,3))+xp(i-1,1); 
s2=(xp(i,2)-xp(i-1,2))/(xp(i,3)-xp(i-1,3)); 
y2(j)=s2*(xpo(3)-xp(i-1,3))+xp(i-1,2); 
plot(y1(j),y2(j),'.') 
hold on 
end 

if xp(i,3)==xpo(3) 

j=j+1; 

yi(j)=xp(i,i); 

y2(j)=xp(i,2); 
plot(y1(j),y2(j),'.') 
hold on 
end 
end 

xlabel('x');ylabel('xdot');title('Problem 14-45:Poincare Section') 
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COA/t~ = 

function xdot=p1445d(t,x) 
global m c k mu F w 
xdot=[x(2); 

(1/m)*(-c*x(2)-k*x(1)-mu*x(1) A 3+F*cos(x(3))); 

w]; 


io~io 
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Problem 14-45:Poincare Section 
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